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Abstract Galaxy clusters, the largest clearly defined objects in our Universe, are ideal lab- 
oratories to study in detail the cosmic evolution of the intergalactic intracluster medium 
(ICM) and the cluster galaxy population. For the ICM, which is heated to X-ray radiating 
temperatures, X-ray spectroscopy is the most important tool to obtain insight into the struc- 
ture and astrophysics of galaxy clusters. The ICM is also the hottest plasma that can be well 
studied under thermal equilibrium conditions. 

In this review we recall the basic principles of the interpretation of X-ray spectra from 
a hot, tenuous plasma and we illustrate the wide range of scientific applications of X-ray 
spectroscopy. The determination of galaxy cluster masses, the most important prerequisite 
for using clusters in cosmological studies, rest crucially on a precise spectroscopic determi- 
nation of the ICM temperature distribution. The study of the thermal structure of the ICM 
provides a very interesting fossil record of the energy release during galaxy formation and 
evolution, giving important constraints on galaxy formation models. The temperature and 
pressure distribution of the ICM gives us important insight into the process of galaxy clus- 
ter merging and the dissipation of the merger energy in form of turbulent motion. Cooling 
cores in the centers of about half of the cluster population are interesting laboratories to 
investigate the interplay between gas cooling, star- and black hole formation and energy 
feedback, which is diagnosed by means of X-ray spectroscopy. The element abundances de- 
duced from X-ray spectra of the ICM provide a cosmic history record of the contribution of 
different supernovae to the nucleosynthesis of heavy elements and their spatial distribution 
partly reflects important transport processes in the ICM. 
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Some discussion of plasma diagnostics for conditions out of thermal equilibrium and an 
outlook on the future prospects of X-ray spectroscopic cluster studies complete our review. 

Keywords X-ray astronomy, Galaxiesxlusters of galaxies, Spectroscopy: X-rays 



1 Introduction 

Galaxy clusters are c ontrary to the ir name more than just a collection of galaxies. In the 
1930s Fritz Zwicky fcwickylll937h discovered that it requires a large amount of unseen 
matter to bind the fast moving galaxies in the Coma galaxy cluster into a long lasting ob- 
ject. Today we have a clear cosmic scenario with galaxy clusters as an integral part of the 
large-scale structure of the Universe. They are the largest matter aggregates within the large- 
scale structure which have collapsed under their own gravity and are closely approaching 
a dynamical equilibrium. Much theoretical effort has been spent to understand and char- 
acterize this equilibrium structure. To the first order, galaxy clusters are now described as 
Dark Matter Halos with a characteristic universal shap e of the Dark Matter potential (e.g. 
iNavarro et al.lll995l : lMoore et alJl999l : lGao et al.l20 08). Observations show indeed, that the 
cluster population can be described as constituting a nearly self-similar family with similar 
shapes of the matter distribution, where the small scatter is due to the different formation 
histories and differently close approaches to the equilibrium configuration. 

Galaxy clusters are therefore very important giant astrophysical laboratories providing 
us with a well characterized physical environment in which we can study many interesting 
astrophysical phenomena and cosmic processes on giant scales (Sarazin 1986). They also al- 
low us to study large coeval galaxy populations and enable us to investigate their evolution in 
connection with the chemical and thermal evolution of the embedding intracluster me dium 
(ICM) (e.g. |Pressledll98d : lDressler et alJll997lJPoggianti et aO 19991 ; iMei et alj|2006h . 

As tracers of the cosmic large-scale structure they are also important probes for cosmol- 
ogy. It is the growth of structure in the matter distribution of the Universe that has a strong 
dependence on the cosmological model parameters and in particular on the nature of Dark 
Matter and Dark Energy. Since galaxy clusters are very sensitive tracers of structure growth, 
a census of the cluster population a s a function of redshift ca n be used to test cosmologi- 
cal models (e.g . Borgani et al.ll2001uSchuecker et al.l l2003a.b; He nry||2004l ; lvikhlinin et alj 
l2003ll2009l : lHenrv et alj2009h . 

The current most important limitation in using galaxy cluster studies for cosmology 
is the calibration of the relation between various observables and cluster mass. Therefore 
a lot of effort is currently being spent to improve the cluster mass determination and the 
under standing of cluster structure (e.g. lAmaudetai]|2007l :l Vikhl inin et ai]|2006l : |PrattetalJ 
2006, 2007). X-ray spectroscopy of the cluster emission plays a crucial technical role in 
this effort to characterize cluster structure precisely, to model the cluster population as a 
family of self-similar objects with explainable deviations, and to establish scaling relations 
of global cluster parameters that allow to draw comprehensive statistical conclusions on 
cluster properties from simple observables. 

Most of the detailed knowledge on galaxy clusters has been obtained in recent years 
through X-ray astronomy. This is due to the fact that the intracluster medium (ICM) has 
been heated to temperatures of tens of Millions of degrees (several keV per particle) which 
causes the hot plasma to emit the bulk of the thermal energy in the regime of soft X-rays. 
Since this is also the photon energy range where the well developed X-ray telescopes come 
into play, galaxy clusters are among the most rewarding study objects for X-ray imaging and 



Fig. 1 The Coma cluster of galaxies as seen in X-rays in the ROSAT All-Sky Survey (underlaying red color) 
and the optically visible galaxy distribution in the Palomar Sky Survey Image (galaxy and stellar images from 
the digitized POSS plate superposed in grey). 

spectro-imaging observations. Fig.Q]shows a composite image of the Coma galaxy cluster, 
where an optical image from the Palomar Sky Survey showing the dense galaxy distribution 
of the Coma cluster is superposed in grey scale on top of an X-ray image from the ROSAT 
All-Sky Survey with X-ray brightness coded in red color. We clearly recognize that the X- 
ray image displays the cluster as one connected entity. This illustrates the fact that galaxy 
clusters are well defined, fundamental building blocks of our Universe. 

As largest cosmic objects they display several interesting astrophysical superlatives: (i) 
the ICM is the hottest thermal equilibrium plasma that we can study in detail, with tempera- 
tures up to two orders of magnitude larger than the temperature in the center of the sun, (ii) 
the gravitational potential of clusters gives rise to the largest effect of light deflection with 
deflection angles e xceeding half an a rcmin, producing the most spectacular gravitational 
lensing effects (e.g. lHattori et al J 19991) . (iii) the hot plasma cloud of the ICM casts the dark- 
est shadows onto the cosmic microwave background through the Sunyaev-Zeldovich effect 
(at wavelengths below about 1.4 mm they are seen as surface brightness enhancements) 
(e.g. lBirkinshawlll999h . and (iv) the merger of galaxy clusters produces the largest energy 
release in the Uni verse after the big bang itself with energies up to orders of 10 63 erg (e.g. 
iFeretti et al.ll2002h . 

X-ray observations and X-ray spectroscopy are the most important tools to obtain de- 
tailed information on cluster properties and the processes occuring in their ICM as will be 
illustrated in this article. In section 2 we will explain how observable X-ray spectra can 
be understood and modeled. In section 3 we show how this spectral modeling is applied 
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Fig. 2 Thermal e quilibrium ionization structure of the oxygen ion family as a function of plasma temperature 
I Bohri ngel 1998ft . 

to study the thermal structure of the ICM that also provides the tool for measuring cluster 
masses. Section 4 deals with the diagnostics of the central regions of clusters where massive 
cooling is prevented by AGN heating in those objects where the cooling time is short enough 
for effective cooling. In section 5 we illustrate how the observed spectral lines can be used 
for the chemical analysis of the ICM and what can be learned from these observations. In 
section 6 we take a look into the physics of plasma under non-equilibrium ionization con- 
ditions and in the last section, 7, we provide an outlook on the capabilities and potential of 
X-ray instruments planned for the future. 



2 X-ray Spectra of Hot Tenuous Plasma 

The plasma of the ICM is very tenuous, with densities of 10~ 5 to 10 _1 cm -3 from the clus- 
ter outskirts to the densest regions of cool core clusters. This low plasma density makes 
the cluster X-ray spectra modeling simple and enables a very straightforward interpretation. 
Three fundamental emission processes involving electronic "transitions" contribute to the 
radiation: free-free or bremsstrahlung radiation caused by the deflection of an electron at 
close fly-by of an ion, free-bound or recombination radiation caused by the capture of an 
electron by an ion following ionization, and bound-bound or deexcitation radiation of an 
electron changing the quantum level in an ion. The first two processes give rise to contin- 
uum radiation and the latter to line radiation. An exception in bound-bound transitions is the 
radiative transition from the 2s to the Is state, which is completely forbidden by angular mo- 
mentum conservation, but can happen as very slow two-photon process. This "two-photon 
radiation" involves a distribution func tion of the branching ratio of the energies of the two 
electrons dSpitzer & Greensteinlll95ll) and thus gives rise to continuum radiation. 

We note that all these radiative processes depend on the collision (or close fly-by) of 
an electron and an ion. Due to the very low density of the plasma all the ions excited by 
collisions have sufficient time for radiative deexcitation before a second deexciting collision 
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Fig. 3 Scheme of the calculation of the radiation spectrum and cooling function of a hot, tenuou s plasma in 
therm al equilibrium. The shaded boxes include the situation of non-equilibrium ionization (from Bbhringer 
119981) . 



occures. Thus contrary to laboratory plasmas, where slow transitions are "forbidden" and 
the corresponding excited states are much more rapidly deexcited by electronic collision, all 
"forbidden" transitions actually happen in the ICM plasma. This leads to a scenario where 
all exciting, recombining, and bremsstrahlung causing collisions lead to the radiation of a 
photon, which is referred to as the thin plasma radiation limit (or "coronal limit", as similar 
conditions prevail in the solar corona). The modeling of the thermal plasma spectrum is 
therefore a book keeping exercise of all the electron ion collisions rates and - in the case 
of several deexcitation channels - their branching ratios. These collision rates are in general 
a function of temperature (for thermal plasma) and the outcome is directly proportional to 
the electron density and the respective ion density. The shape of the resulting spectrum is 
therefore a function of the temperature and chemical composition and its normalization is 
directly proportional to the electron density and the ion density. The latter is true for the low 
plasma density limit that pertains in the ICM. 

We further note that in general all photons so created leave the ICM plasma due to its 
low density. An exception is discussed in section 5.5.3 for the strongest lines in dense cluster 
cores. Thus no radiative transfer calculation is necessary for the interpretation of the X-ray 
spectra of the ICM. This means in particular that the spectrum we observe from a galaxy 
cluster provides an account of the entire ICM plasma, which is very different from e.g. the 
spectra of stars, that provide information on merely a very thin skin on their surface. This is 
another reason why cluster X-ray spectra are so informative and straightforward to interpret. 

As the radiation contribution specific to an ion species is proportional to the density of 
this species in the plasma, the modeling of the overall spectrum also requires the knowledge 
of the chemical composition of the plasma and its ionization structure. For thermal equilib- 
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Fig. 4 Cooling coefficient of hot plasma as a function of the electron temperature from ( Shapiro & Moore 
1976). Two scenarios of the cooling of plasma from an initially hot phase of T ~ 10 7 K are compared. In the 
CIE (collisional ionization equilibrium) case the ionization structure is kept in full thermal equilibrium, in the 
NIE (non-ionization equilibrium) case on the contrary, the ionization structure passively follows the cooling 
plasma and recombination can lag behind the cooling. The latter implies higher ionization stages at the same 
temperature, that is less shell electrons and reduced line cooling. 



rium plasma, the ionization structure is characterized by the balance of the ionization and 
recombination rates for all the ionization steps of an element: 

[I n+ ]n e (Cd + C a d + U" + ] (lH + ]X cel + [He\ He ++ ] X ce2 ) (1) 
= [7<" +1) > e (C rr + C dr ) + [/<" +1)+ ] ([H]Xce3 + [He] X ceA ) . 

Several processes contribute to the ionization and recombination rate. The ionization com- 
prises direct electron impact ionization, C c ,, electron impact excitation into an auto-ionization 
state with subsequent auto-ionization, C ( „, and ionizing charge exchange with H + , He + , 
and He ++ with the rates X ce \ and X ce 2, respectively. Direct ionization is the dominant pro- 
cess. Recombination rates include radiative recombination, C rr , dielectronic recombination, 
Cd r , and electron capture charge exchange with H and #e,with the rates X ce j, and X ce s,. In 
the dielectronic recombination process the recombination collision complex is energetically 
stabilized by the excitation of a second bound electron which subsequently radiates the ex- 
cess energy. This is a resonant process which exhibits a strong energy dependence which 
is most important for ions with many electrons (e.g. Fe ions). The ionization structure is 
calculated from the complete set of linear equations for all important i ons and elements (e.g. 
lArnaud & Rothenfludll985LlArnaud et al.|[l992LlMazzotta et al.lll998h . Fig. |2] shows for ex- 
ample the thermal equilib rium ionization structure of oxygen in the temperature range 10 4 
to 10 8 K l lBohringer|[l998h . 
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Fig. 5 Cooling rate of hot plasma as a function of the plasma te mperature. The contribution to the cooling 
by the ions of different important abundant elements is indicated i Bohringer & Hensler 1989). Most of this 
contribution is in the form of line radiation, which is by far the dominant form of radiation in the temperature 
range from about 10 4 to 2 X 10 6 K. At higher temperatures over most of the regime of interest for the ICM 
the Bremsstrahlung contribution dominates. 



The collision rates that lead to the emission of radiation are then all of the form: 



R = tigiijCit 



nE_ 



n. 



C,,(T), 



(2) 



where he, fii, and tin are the number densities of the elements, ions, and hydrogen nuclei, 
respectively, and Cj x are the relevant collision rate coefficients for ion /. Consequently, \^\ 

represents the fractional abundance of ionization stage /, the relative elemental abun- 
dance, and [^] the Hydro gen nuclei to electron ratio, which is about 1.2. Thus all rates 
are proportional to the squared plasma density, n 2 e . The normalization of the spectrum is 
therefore determined from the emission measure, E, given by: 



dV. 



(3) 



We note that this definition of emission measure is not unique. It conforms to most defi- 
nitions in text books, but the major public plasma radiation codes use the definition E = 
J n e n H dV. The strategy of calculating the radiation spectrum of hot, thin plasma in ther- 
mal equilibrium is thus summarized in Fig. [3] The notion of thermal equilibrium involves 
a thermal Maxwellian equilibrium of the electrons and the thermal ionization equilibrium 
of the ions. A deviation of the temperature of the electrons and the ions will not affect the 
components of the radiation spectrum, but will be detected through thermal line broadening. 
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With a change of the temperature of the plasma the ionization structure is adjusted on 
the time scale of the inverse ionization and recombination rates and thus lags behind the 
temperature change. In certain cases this time lag is significant and has to be included in the 
modeling, as sketched in Fig. [3] by explicitly following the time evolution of the ionization 
structure. As an exa mple, Fig.[4]shows the c alculation of a passively cooling plasma starting 
at high temperature ( Shap iro &Moorel 19761) where the instantaneous ionization equilibrium 
assumption is compared to the full time dependent ionization structure calculations. We 
see that strong deviations occur at temperatures of the order of 10 5 K, but no significant 
deviation from the ionization equilibrium calculations are observed above 10 6 K, which is 
the temperature range in clusters. So far no observed X-ray spectrum of the ICM has required 
a non-ionization equilibrium treatment. Fig.[5]shows the contributions of continuum and line 
emission and the respective contributions of the most important elements (assuming solar 
abundances) to the total radiation power. The regime in which the recombination lags behind 
cooling can easily be identified with the temperature range where line radiation from the 
important, most abundant elements (particularly Fe) is strongly boosting the line cooling. 
We note that continuum bremsstrahlung with a temperature dependence of roughly oc T 1 ^ 2 
is the dominant radiation process at ICM temperatures. 
Early modelin g of the emis s ion spectrum of hot, thi n 



Cox & Tucker 



lasma incl u de: ._ 
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1990h:lLandini & Fossil dl99ll):lGaetz & Salpetej dl983l):lRavmondl dl988h; iMasail Jl984h; 
Masai & Kato ( 1987); Bohrineer & Hensler ( 1989); Schmutzl er & Tscharnuterl l ll993l) ; ISutherland & Dopital 
{ 1993h : Brickhouse et al. ( 1995S): iRavmond & B rickhousd dl996h. Now the two most fre- 
quently used plasma radiation codes are the MEKAL code (iGronenschild & Mewejjl97ij ; 
Mewe & Gronenschilj|l98ll : lMewe et al.ll 1 985l : iMewdl 1 99dl . 1 199 ll : iKaastra & MeweMl994l ; 
Mew e et al.ll 19951) and the APEC code (Smith et al. 2001). While the older literature con- 



tains all the necessary information about the radiation codes including all the details of the 
rate calculations and much detail on the underlying physical processes, the modern codes 
are too complex to be easily described in the literature. Thus, unfortunately very little doc- 
umentation about the new development of the radiation codes is available in written form 
after the mid 1990s. 

Fig.[6]shows typical spectra at temperatures of 10 5 , 10 6 , 10 7 , and 10 8 K with solar ele- 
ment abundances, indicating the contribution by line radiation and the continuum emission 
from bremsstrahlung, recombination and two-photon transitions. We clearly see the increas- 
ing dominance of bremsstrahlung with increasing temperature, which reflects the fact that 
fewer ions retain electrons and the plasma is almost completely ionized at the higher tem- 
peratures. 



3 The Study of the thermal structure of the intracluster medium 

We have seen in the previous section that the shape of the spectrum for a thermal equilibrium 
plasma is determined by the plasma temperature and the elemental abundances. This is 
therefore the basic information we derive from the spectral analysis of the ICM radiation: 
a temperature measurement and a chemical analysis. We consequently illustrate in this, and 
the next chapter, the scientific insights gained from temperature measurements from the 
state-of-the-art spectral analysis, and in section 5 the lessons learned from the chemical 
analysis of the ICM. 
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Fig. 6 X-ray spectra for solar abundance at different plasma temperatures. The continuum contributions from 
bremsstrahlung (blue), recombination radiation, characterized by the sharp ionization edges (green), and 2- 
photon radiation (red) are indicated. At the highest temperatures relevant for massive clusters of g alaxies 
bremsstrahlung is the dominant radiation process (from the work described in Bohringer & Hensler 1989). 
The major emission lines in the panels for the higher temperatures relevant for galaxy clusters are designated 
by the elements from which they originate (The labels Fe-L ans Si-L refer to transitions into the L-shell in 
ions of Fe and Si, respectively, and two other lines with roman numbers carry the designation of the ions from 
which they originate involving transitions within the L-shells. 



It is mostly the shape of the continuum spectrum, dominated by bremsstrahlung, that 

provides information on the temperature. The spectral energy distribution for th e thermal 

brem sstrahlung spectrum for the collision of an electron with ion, is given by (e.g. lGronenschild & Mewel 
Il978h : 

16 e 6 ( In \ 1/2 , / -hv \ 

e(v) = j — n e m Z 2 g ff {Z, T x , v) exp — — , (4) 

3 m e c z \ 3m e k B T x J \k B T x / 

where m e and n e are the electron mass and density, respectively, n ; is the respective ion 
density, Z is the effective charge of the ion, and gff is the gaunt factor, a quantity close to 
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Fig. 7 X-ray spectrum of the ICM plasma in the halo of M87 at a temperature around 1 .7 keV observed 
with XMM-NewtonEMOS, showing exclusively the spectral range with the H- and He-like lines of Si and S 
I Matsushita et al. 2002). The spectrum is fitted with a one-temperature MEKAL model (black line) and a two- 
temperature model with the temperature components 1.4 keV and 3.3 keV (grey line). The one-temperature 
model provides a better fit, supporting the view that the temperature distribution is very narrow. For more 
details see the original publication. 



unity which must be calculated numerically through a quantum-mechanical treatment. The 
most prominent spectral signature of Tx is the sharp cut-off of the spectrum at high energies, 
due to the exponential term with the argument -hv/k B T. As long as this cut-off is seen in the 
energy window of the telescope, one has a good handle on the temperature measurement. 
The element abundances are mostly reflected in the intensity of the spectral lines. Only 
in well resolved spectra of very high quality, where we can observe several lines of the 
same element, is the temperature also constrained by line strength ratios. A prominent such 
case occurs at temperatures around 2-3 keV where we can observe the K-shell and L-shell 
lines of iron simultaneously in high photon count spectra. Another example of temperature 
diagnostics based on the study of the hyd rogen and helium like Si and S lines is shown in 
Fig.Q(taken from lMatsushita et alj|2002l) . This work provides a further nice illustrations of 
how different temperature indicators can be used to check the consistency of the temperature 
measurement and to test if the plasma has a distribution of temperature phases. 

The first X-ray spectral observations of clusters lacked sufficient angular resolution and 
thus involved the total ICM emission of the targets and only provided us with information 
on the bulk temperature of the cluster ICM. What does this parameter tell us about a clus- 
ter? Galaxy clusters form from the gravitational collapse of overdense regions in the matter 
density distribution in the Universe and subsequently approach an equilibrium configuration 
characterized by a virial relation: 
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Fig. 8 Correlation of the ICM temperature with the velocity dispersion of the galaxies in clusters I Mushotzkv 
1984). The temperatures have been derived from X-ray spectra taken with the HEAO-1 satellite A-2 experi- 
ment. The dots are data for non-cD clusters and the open symbols for cD clusters (clusters with a centr al dom- 
inant, cD type galaxy). The lines show predictions of polytropic models with various indices (see Mushotzkv 
[19841). 



where the mass, M, refers to the total mass of the galaxy clusters including the dark matter. 
Analogously to the virial equilibrium of galaxies and dark matter particles, the ICM plasma 
thermalizes and attains a "virial temperature" which reflects the depth of the gravitational 
potential of the cluster. In the collapse process the potential energy of the ICM is converted 
to internal heat. If the gravitational potentials of clusters of different mass hav e a self- similar 
shape , as implied by nu merical simulations of gravitational collapse (e.g. iNavarro et alj 
ll995l : lMooreetal.lll999h . then one finds the following self-similar relation between clus- 
ter mass and ICM temperature: 

T cc cr 2 DM oc M/R cc M 2/ \ (6) 

where cr DM is the velocity dispersion of the dark matter particles. (The rightmost relation is 
due to M oc R 3 and pom = const, in the self-similar model). 

The analysis of the first galaxy cluster spectra indirectly confirmed this trend by show- 
ing that the ICM temperature increases with th e galaxy velocity dispersion. Fig. [8] shows 
one of the first such relations dMushotzky|ll984) with X-ray spectral observations from the 
HEAO-1 A2 satellite experiment. Modern versions of these relations - given here directly 
as M — Tx relations - are shown in Fig. [9] f or selected r egular clusters from lArnaud et al.l 
d2007l) and for cluster sample from iKotov & Vikhlinir] d2006h . respectively. The spectro- 
scopically determined ICM temp erature thus turns out to be one of the best mass proxies as 
single observable parameter (e.g. iKravtsovet al] |2006l) . The tightest relations are obtained, 
if the core regions are excluded in the global temperature measurement, due to the dispro- 
portionate influence of the central cool cores as will become aparant below, and it has thus 
become standard to quote the mean temperature in the radial region r = 0.15-lx 7-500 as the 
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Fig. 9 Left: The mass - temperature relation of a sample of regular clusters from I Arnaud et al. 2007}. The 
best fitting slope (solid line) is 1.71 ± 0.09 slightly steeper than the self similar relation. Right : The mass 
- temperature relation of galaxy clusters and its evolution with redshift derived from two cluster samples at 
< z >~ 0.05 and < z >~ 0.55 iKotov & Vikhlinin 2006). The relations have been displaced by a factor of 
10 for better visibility. The model predicted redshift evolution (clusters of equal mass get hotter with look- 
back time) has been accounted for by scaling with the factor h(z). This predicted evolutionary trend is well 
supported by the data. 



most reliable single observable mass proxy for clusters (e.g. lAmaud et alj|2005l : [Pratt et alj 

boosts . 

The advanced X-ray observatories Chandra and XMM-Newton now routinely also pro- 
vide localized measurements of the ICM temperature. First attempts had already been made 
using observations with the ROSAT o bservatory e.g. to produce temperature map s featuring 
still large temperature uncertainties dBriel & Henrvl U994k iHenrv & Briellll995t) and later 
by means of the ASCA observatory, but due to its comparably low angular resolution most 
of th e results remained somewhat ambiguous and depende d on the analysis method used 
(e.g. Ilkebe et al.l fl997l: iMarkevitch et alJl 19981 : IWhitd boOO). With some studies by means 
of the Italian-Dutch X-ray mission BeppoSAX some overview on the temperature structure 
of clusters was obtained, with decreasing temperature profiles outside ~ 0.2rsoo and a di- 
versity of the behaviour in th e core region depending on the thermal history of the cluster 
( iDe Grandi & Molendill2002h . This picture has become very precise with systematic studies 
with XMM-Newton and Chandra. Fig. fTOlshows results on the temperature profil es of cluster 
samples studied with both observatories dVikhlinin et ap2006l ; IPratt et al]|2007l) . 

To better understand these results, we have to make a brief excursion to scaling relations. 
For the comparison of the structure of galaxy clusters of different mass based on the model 
of self-similarity mentioned above, we need a fiducial radius which identifies corresponding 
scales in clusters of different mass. A useful definition of such a scale is the radius at which 
the mean mass density of the cluster is larger by a certain factor than the critical density 
of the Universe at the redshift of the cluster. One of the justifications of this picture comes 
from the spherical collapse model in a critical density universe. A good practical choice for 

1 r5oo is the radius where the mean total density of the cluster is 500 times the mass density of a Universe 
with critical density. This radius describes the same relative scale of clusters of different mass in the self- 
similar scenario. See also explanation later in the text. 
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Fig. 10 Scaled temper ature profiles of gala xy clusters derived from spectroscopic observations with the 
Chandra satellite (left, Vikhlinin et al. 2006) and results obtained for a sample of clusters observed with 
XMM-Newton (right. lPratt et alj2007l) . 



the overdensity factor is 500 (see e.g. Evrard et al 1996 who show that inside a region of 
~ r 5 ()o randomized galaxy and mass particle orbits dominate clearly over infalling material), 
not least from an X-ray observers point of view, since now a larger number of high quality 
Cha ndra and XMM-Newton data cover the c luster ICM spectroscopically out to this radius 
(e.g. IVikhlinin et all EoOfil: IPratt et all Eobl). We will most often make use of this scaling 
radius of r^oo, but also use other values of overdensity depending on the literature from 
which we draw the examples. 

In Fig. [TO] the temperature profiles shown are scaled by the overdensity radii (with dif- 
ferent overdensity values, which has no importance here). We note that there is a large 
diversity of the profiles at small radii. There are two classes of galaxy clusters: the clus- 
ters with dense ICM cores show temperature profiles with T x decreasing towards the center 
in the core regions, while in clusters with moderate central densities (typically below 10~ 2 
cm -3 ) the temperature profiles are found to be flat or even slightly increasing towards the 
cluster center. The clusters will be called cool core and non-cool core clusters in the fol- 
lowing. Due to the logarithmic display of the radial range, the inner regions appear more 
prominent in Fig. [TTJ Most of the cluster size and volume is in the outer region, where the 
temperature profiles become very similar and appear squeezed into a narrow zone in scaled 
radius. This shows that on large scale the thermal structure of clusters can be viewed as fol- 
lowing closely a self-similar model. Characterizing this model and understanding the scatter 
around the mean model is one of the current important observational goals of X-ray cluster 
research. 

One of the difficulties in deriving the temperature distribution of the ICM in clusters 
from X-ray spectroscopic observations is the fact that the observed radiation is the result 
of an integral of radiative emission along the entire line-of-sight through the cluster. There- 
fore one observational task is to "deproject" the cluster spectra, which can only be done 
by assuming a certain three-dimensional geometrical shape of the cluster, in general spher- 
ical symmetry. The requirements for successful spectral deprojection are very good photon 
statistics and an angular resolution of the data (instrument point spread function) that is bet- 
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ter than the corresponding radial binning used. In general the deprojected spectra are much 
more noisy than the observed spectra, but most often the deprojection results show only 
small specific differences to the unprojected data (because most of the emission shaping the 
spectrum comes from the center-most region). 

A further problem in deriving the temperature structure from spectroscopy arises if the 
plasma has more than one temperature phase locally, that is within the observed patch from 
which the spectrum is extracted and within the deprojection bin. The question is then: what 
determines the resulting temperature when this composite spectrum is fit with a single-phase 
spectral model? A first intuitive guess is to assume that it is an average of the temperatures 
weighted by the radiative emission co ntributions of the different p hases, the so-called emis- 
sion weighted t emper ature, T ew (e.g. iMathiesen & Evr ard 2001). However, as shown by 



iMazzotta et al.1 (2004) the resulting fitted temperature is generally biased low compared to 
T ew . Mazzotta et al.l 1 2004T) derive a simple analytical approximation to calculate the fitting 



result from a temperature mixture with an accuracy of a few percent for ICM temperatures 
above 3 keV. The resulting fitted temperature they call spectroscopic-like temperature, T s {. 
This approach is based on the fact that the temperature dependence of the spectra according 
to Equ. 4 can be expressed as: 

e(E, T) ~ n] &T,m) T~ 1 ' 2 «p(~ J;) ~ t(T,m) T~ 112 (l - ^ (7) 

where E is the photon energy, m is the metallicity (heavy element abudance) of the 
plasma, and f (T, m) includes the temperature dependence of the gaunt factor and some effect 
of emission lines. The right hand side is a Taylor series expansion for low enough energies 
as covered by the energy window of X-ray telescopes (for k B T > 3keV). Then it is easy to 
show that temperature phase addition with subsequent single temperature fitting leads to the 
following expression for the spectroscopic like temperature: 

r n 2 T "/ T i/2 dv 

" ~ Jn 2 T"/Ty 2 dV 

with tests showing that a = 0.75 is a good approximation. IVikhlinirj J20061) has gen- 
eralized this approach also to low temperatures using a method that requires some more 
complex numerical simulations that take into account the specifics of the spectroscopic in- 
strument, energy band and galactic absorption providing satisfying results with errors of a 
few percent. 

The latter approach allows us to derive T s [ for a known temperature distribution, but 
provides no means to correct for a temperature bias, if the distribution is unknown. For 
example the determination of cluster masses (as explained below) rests on a local average of 
the pressure which requires a mass weighted temperature average for known density. Since 
T;/ is always biased low with respect to mass averaging, the mass will be unde restimated in 
the presence of u nresolved multi-temperature structure (e.g. lRasia et "aT] |2005h . A study by 
Rasia etal. (2006) showed that the temperature bias can easily lead to a mass underestimate 
by 10% on average. With current instrumentation it is difficult to unveil multi-temperature 
structure in the temperature range above 3 keV. A simple exercise that can be performed with 
XSPEC is to take an equal emission measure of 4 keV and 8 keV plasma and to simulate an 
analysis with the instruments of XMM-Newton. Even with a simulated exposure that gives 
several ten million photons, a single temperature model with T ~ 5.5 keV gives a perfect fit 
and we are unable to recover the original input phases. At lower temperatures and very good 
photon statistics more diagnostics can be performed as seen e.g. in the study shown in Fig. [7] 
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Fig. 11 Upper Left: Mass profil es for a sample of regular c lusters with a range of temperatures derived from 
observations with XMM-Newton ( Pointecouteau et al. 2005). Upper Right: Same mass profiles as in the left 
panel scaled in radius by rim and in mass by M200, the mass at the radius of r2oo- The solid line shows the 
fit of a NFW mass distribution model which provides a good description of the data. Lower Panel: Scaled 

[ as density profiles an d scaled total density profiles for a sample of galaxy clusters observed with Chandra 
Vikhlinin et al. 2006). The yellow line shows the best fitting NFW-model for the total density profile. 



I Matsushita et al. 2002 ) and simila r studies in other cool core regions e.g. by iFabian et alj 
(2005) or lSimionescu et al. (2008a) shown below, due to the emisson lines characteristic to 
the different temperature phases. 



3.1 Cluster Mass Determination 



One of the most important applications of the knowledge of the temperature structure in 
the ICM is the mass determination of galaxy clusters. Assuming the ICM is in hydrostatic 
equilibrium and the cluster has approximately spherical symmetry, the distribution of the 
total mass in the cluster is given by the density and temperature profiles by the following 
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equation: 

M W = -^(^ + ^£), (9) 
ljm p \ logr logr ) 

where G and k B are the gravitational and Boltzmann's constant, respectively, /i is the mean 
particle mass (~ 0.6), m p the mass of the proton, and p g the gas density. 

One of the generally interesting results from recent efforts of precise mass determina- 
tion is the conclusion that the "NFW" dark matter halo model proposed by Navarro, Frenk 
& White (1995,1997) gives a consistent description of the mass distribution of seemingly 
relaxed clusters (clusters that show a high degree o f symmetry in projection and no obv ious 



signs of recent merger activ i ty) as shown in Fig 
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(2006)). First attempts have been made of a rigorous comparison of the cluster masses de- 
term ined from an X- ray analysis and from gravitational lensing studies dZhang et al.ll2007l : 
Mah davi et al.ll2008h showing that the X-ray mass determination provides a much smaller 
individual uncertainty. We also have first indications that the lensing masses are generally 
higher by about 12(±15)% which may partly reflect the above mentioned temperature bias 
and additional unaccounted turbulent pressure of the ICM if we assume that the lensing 
masses are essentially unbiased. This mass calibration by combining several methods will 
constitute an important effort in the coming years. 



3.2 Fossil Record of Structure and Galaxy Formation in the ICM Thermal Structure 

The temperature distribution in the ICM is also the key to the characterization of the thermal 
history and thermal structure, most conveniently expressed by the entropy structure of the 
cluster's ICM. The definition of "entropy", S , as used in the astrophysical literature of galaxy 
clusters deviates from the general definition in physics: 

S = k -%- do) 

n 

Th is definition is related to the general definition of entropy, s, by s = k B In (s^^Qxrrip) 5 ^ + 
so dVoitll2005h . The S value can be seen as a parameter to lable adiabates, and thus S stays 
fixed in any hydrodynamic evolution of the ICM in which all processes are adiabatic; S 
increases e.g. in shock waves or with dissipation of turbulent motion and descreases with 
radiative cooling. 

The most important application of studying the entropy structure in galaxy clusters is the 
distinction between heating by gravitational and non-gravitational processes. The term grav- 
itational processes refers here to heat that comes from conversion of potential to heat energy 
in the formation of structure and the collapse of clusters. This form of heat is determined 
by structure formation processes that are mainly driven by the dark matter in the Universe. 
Non-gravitational processes are connected to energy input into the ICM in the form of star 
formation driven galactic winds (which at the same time lead to chemical enrichment of the 
ICM as discussed in section 5) and by AGN providing energy input via jets of relativistic 
plasma. 

The way to distinguish the two forms of heating is the study of scaling relations of 
thermal cluster properties. The gravitational heating of the ICM during cluster formation 
is proportional to the depth of the gravitational potential well of the cluster and thus the 
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Log 10 [r(keV)] 

Fig. 12 Entropy in grou ps and cluste rs measured at a fiducial radius of 0. 1 r v , rla / as a function of the system's 
ICM temperature (Ponman et al. 1999). It was inferred from this plot that the entropy appears to converge at 
low Tx to a characteristic floor value of ~ 100 h ' keV cm 2 . 



entropy gain is proportional to the ICM temperature. Heating processes connected to feed- 
back from the galaxy population are expected to be related (proportional) to the total stellar 
mass in the cluster. Assuming to first order a constant stellar to gas mass ratio of the bary- 
onic cluster component, the amount of non-gravitational energy input per unit gas mass is 
than also to first order constant. Therefore in this most simple consideration the gravita- 
tional entropy contribution is proportional to the virial and ICM temperature, while the non- 
gravitational entropy contribution is an additive const ant. The evidence for non-gravitational 
entropy characte rized in this way w as provided by IPonman et al. I dl999l) as illustrated in 
Fig. \M (see also lDavid et~ai1 fl996) who showed with ROSAT observations of cooler sys- 
tems that groups have flatter entropy profiles than cluste r s). The results imply a so -called 
entropy floor of about 135 keV cm 2 jPonman et alJI 19991 : iLlovd-Davies et al.ll200(t) . More 
modern results do not show the sudden appearance of a step, but rather a slope of the entropy 
- temperature relation which is shallower than the slope of 1 expected from purely gravita- 
tional models: instead of the gravitational slope of S{r scaled ) oc T, relations around oc T 1 ^ 
are fou nd. Purely gravitational models of smooth and cold accretion of gas into the clus- 
ter (e.g. lTozzi & Normanll200lh predict entropy profiles in clusters of the form S{r) oc r 11 . 
Observations show profiles which can be as steep as this prediction for the most massive 
systems, but they are shallower for smaller clusters, which is again a manifestation of the 
fact that non-gravitational processes have a larger effect on systems with a smaller mass 
and thus a s hallower gravitational potential. As a representative result we show the study of 
IPrattl d2006l) involving 10 relaxed appearing clusters (Fig. 113b. which illustrates the radial 
and ICM temperature scaling behaviour of the cluster ICM entropy. The d iagnostics of en- 
tropy profiles has now been pursued in detail for larger cluster samples by Cavagnolo et alj 
(2009) and representative cluster samples by Prat et al. (2009, in preparation) providing a 
more detailed picture, the proper interpretation of which has to be assisted by simulations. 
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Fig. 13 Left: The entropy - temperature relation for 10 regular galaxy clusters determined at four different 
fiducial radii I Pratt 2006). The slopes of the relations as determined from the orthogonal BCES method 
iAkritas & Bershadvll 19931 are for increasing radii: 0.49 ± 0.15, 0.62 ± 0.11, 0.64 ± 0.11, and 0.62 ± 0.08. 
Right: Scaled radial entropy profiles for the same 10 clusters. The radius is scaled by noo- The entropy, S , has 
been scaled with an entropy-temperature relation proportional to T° 65 using the global cluster temperature. 
The grey shaded are a corresponds to the lcr standard deviation. The dashed line shows radial dependence 
S oc r 108 l:Pratt 20061). 



Theoretical modeling to explain the quantitative entropy scaling behaviour involves 
"preheating" (early entropy increase of the intergalactic medium due to star formation be- 
fore the formation of the cluster), cooling and condensing of the low entropy material which 
increases S in the remaining gas phase, and feedback heating processes in the cluster. A 
complicated balance of all these processes seems to be necessary to reproduce the entropy 
profiles, the scaling relations, and the amount of baryons that are converted into stars in the 
cluster volum e and efforts to obtain a completely consistent and satisfying model are sti ll 
ongoing (e.g. lBorgani et al.ll2004l200l ; IVoit et alJl2003l ; l\^l2005l ; lMcCarthv et al.ll2008h . 



3.3 At and beyond the virial radius 

The low, stable, and well understood instrumental and particle background properties of the 
Suzaku satellite all owed to perform the fi rst pioneering studies of clusters out to and beyond 
their virial radius. iReiprich et all (2008) determined the temperature profile for the rich, 
massive, cooling core cluster Abell 2204 from ~ 10-1 800 kpc, close to the estimated r2oo- 
They find that the temperature profile between 0.3-1.0^200 is consistent with a drop of 0.6, 
as predicted by simulations. The first temperature, density, and entropy profile s beyond ^200 
were reported for another rich, massive cooling core cluster PKS 0745-191 ( George et al. 
2009). Between 0.3-lr2oo they measure an average temperature drop of ~ 70%. They find 
that near to the virial radius the observed entropy profile is lower than that expected for the 
heating by gravitational collapse. Fig.[l4]shows the observed profiles in four directions and 
the average profiles. The dotted curve at the bottom shows the expected entropy profile for 
heating by gravitational collapse, S °c r , and the vertical dashed line shows the estimated 
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Fig. 14 Temperature, density and entropy profiles of PKS 0745-191 determined by Suzaku out to and beyond 
the virial radius in each direction. The dotted curve at the bottom shows the expected entrop y profile for 
heatin g by gravitational collapse S oc r . The vertical dashed line shows the estimate of r2oo I Geor geet all 
120091) . 



7"200- The fall of entropy beyond raoo in the NW direction is interpreted by the authors as 
evidence for an accretion shock from cooler material falling on the cluster along a filament. 

X-ray emission from hot plasma associated with a filament connecting two massive 
clusters of galaxies , Abell 222 and Abell 223, was recently detected using XMM-Newton 
dWemeretal1 l2008). The detection of the tenuous gas permeating the filament was pos- 
sible because of its favorable orientation approximately along our line-of-sight. The tem- 
perature of the detected gas is kT = 0.91 + 0.25 keV and assuming that the length of the 
filament along the line-of-sight is / = 15 Mpc its baryon density is ~ 3.4x 10~ 5 cm -3 , which 
corresponds to a baryon over-density of * 150. The entropy of the gas in the filament is 
S m 870 keV cm 2 Z 1 ^ 3 , which suggests strong preheating. We note that if the detected emis- 
sion would be associated with the outer region of radially asymmetric clusters (deriving 
the gas density by means of the usual geometric deprojection), the corresponding entropy 
of ~ 420 keV cm 2 would be much lower than that expected for the virialised cluster gas 
(~ 1000 keV cm 2 ) which has already passed through an accretion shock. 
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Fig. 15 Temperature map of the merging cluster A3921 by (Belsole et al. 2005). The temperature map has 
been produced from XMM-Newton observation s obtained with the detectors EMOS1 and EMOS2 with the 
multi-scale spectro-imaging technique of Bourdin et al. 12004). The barlike hot region to the NE is interpre- 
tated as the signature of an off-axis merger a t a stage short after the first close encounter similar to features 
observed in simulations iBelsole et al. 2005; Ricker & Sarazin 2001). Superposed to the temperature map 
are the contour lines of the surface brightness distribution in the 0.3 - 10 keV band which was adaptively 
smoothed. The contours are spaced logarithmicly. 



3.4 Diagnostics of Cluster Mergers 

Galaxy clusters grow throughout the recent history of our Universe by accreting matter 
from their surroundings, preferentially by clumpy accretion from the cosmic network matter 
filaments in the intersections of which they are embedded (e.g. iBraglia et al. I l2007h . From 
time to time major merger events happen in this accretion process, where two larger systems 
are attracted to each other and merg e violently. Thes e cluster mergers have always been very 
attractive study objects (see e.g. iFeretti et al.ll2002l . for a review). The thermal structure of 
the ICM again offers good perspectives to unveil the merger configuration and its history as 
well as to understand the heating processes of the ICM resulting from the energy release of 
the merger. 

We will illustrate the diagnostic potential of X-ray spectroscopy in these studies by con- 
centrating here mostly on the example of the most dramatic merging system, 1E0657-56. 
But before we describe an earlier result obtained with XMM- Newton for the diagn ostics 
of the off-axis merger in the cluster A3921 shown in Fig.[T5lbv lBelsole et al.1 d2005h . The 
temperature map has been obtained by ap plying the mul t i-scale spectro-imaging technique 
based on a wavelet analysis described bv iBourdin et alj J2004I) . The general result of this 
analysis was subsequently confirmed by spectroscopy of the regions highlighted in the tem- 
perature map. While the general gas temperature of the undisturbed Eastern region of the 
main cluster has a temperature around 4.9 keV, a hot, bar like region with an extent of about 
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Fig. 16 Thermodynamic maps for the ICM of the "bullet cluster", 1E0657-56 jMillion & AllerJl2008h pro- 
duced from a 500 ksec observation of this cluster with Chandra. The panels show the X-ray surface brightness 
in the 0.8 to 7 keV band (upper left), the temperature, kgT, in units of keV (upper right), the projected pres- 
sure in units of keV cm~ 5 ' 2 arcsec -1 (lower left) and projected entropy in units of keV cm 5 ' 3 arcsec -1 (lower 
right). The shock front preceding the bullet is marked by an arrow in the temperature map. 



160 x 280 kpc is observed in the highly distorted North- Western side which has been found 
to have a temperature aroun d 7.25 keV. The signatures have been interpreted as an off-axis 
merger dBelsole et alJ2005h . The observed features have been found to be similar to those in 
the simulations of an off axis merger with a mass ratio of 1 :3 an impact parameter of b = 5r s 
(where r s is the scale radius of the NFW profile describing the density distribution of the 
cluster) a nd a time of less than ab out 1 Gyr after the X-ray luminosity reaches its peak as 
shown bv lRicker & Sarazinl l200lh in their Fig. 7. This system seems to show the typical ap- 
peara nce of an intermediate stage off-axis merger similar to other cases (e.g. IReiprich et alj 
l2004h . 

The most detailed and interesting data set on a merging cluster system is that of 1E0657- 
56 at a redshift of z = 0.297. This cluste r was observed in a very deep observation with 
Chandra with an exposure of 500 ksec I IMarkevitchl 12006) and also with XMM-Newton 
(Zhang et al. 2006, Finoguenov et al. 2005). As a first exercise to get an overview on the 
thermodynamic structure of the cluster ICM, it has become standard for such deep, high 
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Fig. 17 Left: Surface brightness profile across the shock front preceding the bullet I Markevitch et al. 2002). 
The solid line shows the best fit model for the three-dimensional density distribution with a sharp spherical 
discontinuity at the shock. Right : Projected temperature profile in a narrow sector across the shock front. 
The vertical l ines show the bou ndary of the bullet and the shock. The dashed line shows the average pre-shock 
temperature ( Markevitch 2006). Both data sets have been derived from deep Chandra observations. 



photon statistic observations to produce maps of the temperature and density distributions 
of the ICM and to infer the pre ssure and entropy distri bution from these quantities. Fig. [l6l 
shows such maps produced by ( [Million & AllerfeOOSt) . The temperature is determined from 
spectroscopy of about 100 image pixels, each containing at least 3000 p hotons. The bi nning 
of these pixels has been guided by regions of similar surface brightness ( Sanders 2006). The 
ICM density is derived from the surface brightness distribution by assuming a certain geom- 
etry to deproject the emission distribution along the line-of-sight. Although the so derived 
temperature in the temperature map is a projected quantity, it still gives a good impression 
about the temperature in the central bin in the line-of-sight because of the large weight of 
the innermost bin due to the square density dependence of the emissivity and the steep den- 
sity profile. From these maps of approximate temperature and density distribution in a cross 
section of the cluster, the approximate distribution of pressure, P = nk B T x , and entropy, 
S = T/n 2 / 3 , can be constructed. The surface brightness image in Fig.[l6]shows a disturbed 
larger cluster component and a compact, cone like structure to the West. The latteris iden- 
tified with a compact subcluster, flying through the main cluster at high relative velocity, 
for which optical images provide further evidence. This component has been named "the 
bullet" from which the popular name of the whole system, "the bullet cluster", originates. 
In the temperature map the bullet shows up as a cool core which also has the lowest entropy 
in the entropy map. The Mach cone like shape of the bullet suggest that it may fly with 
supersonic velocity. In the temperature map we observe that the region in front of the bullet 
has a strongly enhanced temperature that goes hand in hand with high pressure and elevated 
entropy. This is the signature of a region heated by a shock preceding and being detached 
from the bullet. The high temperature/entropy region shows a sharp edge at the shock. In 
the pressure map the region of the bullet shows little enhancement (apart from the narrow 
band of shock compressed ICM), while the highest pressure values are found in the center of 
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the overall system structure. Thus, despite the disturbances, the pressure maximum is most 
probably still indicating the region of the deepest gravitational potential. 

To obtain more quantitativ e information about the nature of the shock in front of the 
bullet, dMarkevitch et al.ll2002h has studied the temperature and density distribution around 
the region of the shocked ICM in more detail. Fig.[l7]shows two sharp edges in the surface 
brightness profile in the left panel. The inner bump is the contact discontinuity that separates 
the bullet from the shock heated ICM, the second bump at 90" is the shock. From a density 
model that fits the projected surface brightness one infers a density jump that corresponds 
to Mach number 3 ± 0.4 shock. The corresponding temperature jump as derived from the 
projected spectra, is shown in the right panel. The relatively high Mach number implies 
a relative velocity of the bullet and the ICM of about 4700 km/s. For an explanation of 
this high velocity see recent simulations that try to reconstruct the merger configuration of 
1E0657-56 and recover such high v elocities which are partly boosted by the accretion inflow 
of ICM at the position of the bullet dSpringel & Farrarll2007l : lMastropietro & Burkertll20 08). 
In section 6.2 we discuss the use of the observed shock structure in this cluster to study the 
processes controlling the thermalization of the plasma behind the shock. 

So far only two more cases of clear shock signatures in mergin g clusters are known: 
A 520 dMarkevitch et al.l2005l : lMarkevitch & Vikhlininl2007l) : A22 19 {Million & Allenl2008h . 

Another type of interesting structures in the ICM are the so called cold fronts, which 
also show up in the surface brightness images of clusters as sharp surface brightness discon- 
tinuities. But thanks to the spectroscopically determined temperatures one can show that the 
pressure across the cold fronts is continuous, such that they are boundaries between c older 
denser plasma and a more tenuous, hotter environment. Markevitch & Vikhlinin (2007) pro- 
vide a nice recent review about the X-ray observations and the physics of cold fronts and 
ICM shocks. 



3.5 Observational Studies of Turbulence of the ICM 

Galax y cluster mergers are expected to induce turbulent motions into the ICM dSunvaev et alj 
12003b . An investigation of the presence of stochastic turbulence has been performed by 
m eans of a spectral reduct ion of deep XMM-Newton observations of the Coma galaxy cluster 
bv lSchuecker et alj ||2004f). The turbulence fluctuation spectrum was studied in the projected 
pressure map of the ICM in the central region of the Coma cluster. The Coma cluster has a 
very flat appearan ce, characterized by a very large core radius of the X-ray surface bright- 
ness of ~ 400 kpc I IBriel et all 1992). This enables us to treat the configuration of the central 
region of Coma in the first approximation as a slab geometry, with corrections to the power 
spectrum applied later. The analysis of the fluctuation spectrum testing for a turbulent power 
law spectrum was performed, with the pressure rather than with density or temperature fluc- 
tuations, not to confuse the turbulent fluctuations with static entropy fluctuations in pressure 
equilibrium, which would be characterized by contact discontinuities. In fact, Fig. 1 191 (left) 
shows that turbulent pressure fluctuations dominate also the density fluctuations rather than 
contact discontinuities. 

The map of the projected pressure distribution in the center of the Coma ICM as shown 
in Fig. [T8] was obtained by calculating the gas density from the X-ray surface brightness 
(with an assumed depth of the ICM in the line of sight) and deriving the temperature by a 
spectral analysis of the data in pixels of 20 by 20 arcsec 2 and alternatively 40 by 40 arcsec 2 , 
yielding the pressure by means of the ideal gas equation of state. We see the pressure fluctu- 
ations clearly in Fig. [18] A Fourier analysis of these fluctuations seen in projection results in 
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Fig. 18 Map of the projected pressure distribution of the ICM in the central region of the Coma galaxy 
cluster. The scale size of 14 5 kpc corresponds to th e largest turbulent eddies traced by the observed pressure 
fluctuation power spectrum ( Schuecker et al. 200J). 



the power spectrum shown in Fig. [19] (right). This spectrum has been corrected for the con- 
tribution of Poisson noise and the overall shape of the ICM surface brightness in Coma. The 
observed power spectrum is characterized by a shape very close to a power law. In Fig. [19] 
we show how different 3-dimensional power law functions compare to the observed spec- 
trum if they are projected in the same way as the observed spectrum. The observations lay 
in between a power law exponent of 5/3 and 7/3. An exponent o f 7/3 is the one expe cted 
for the pre s sure fl uctuation spectrum of the classical prediction of Kolmog orovl dl94lh and 
lOboukhovl dl94ll) . The original work by Kolmogorov was developed for an incompressible 
fluid, and in the ICM also magnetic fields play a not fully quantified role, and therefore 
the the case of the cluster ICM is not easily comparable to the classical picture. Neverthe- 
less, a very similar dimensional consideration of turbulence (see e.g. Landau & Lifshitz Vol. 
VI) is surely applicable to the cluster ICM where most importantly a scale free spectrum 
is expected between the driving and dissipation scale. The observed signature of a nearly 
scale free power law power spectrum is thus a very interesting result. The observations have 
therefore been interpreted in the way that a nearly classical turbu lence configuration has 
been established in the Coma cluster ICM dSchueckeretal .112004]) . This can most proba- 
bly be explained b y the fact that Coma is generally believed to be a post-merger cluster 
dWhite et al.lll993h . 

Important information on the amount of turbulence and turbulent pressure support can 
also be obtained by measuring the level of resonant scattering in emission lines observed in 
the ICM. The ICM is generally assumed to be optically thin. This is certainly true for most 
of the emitted X-ray photons, but for the strongest resonance lines the ICM can be moder- 
ately optically thick dGilfanov et aiT l987). What happens is that for strong resonance lines 
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Fig. 19 Left: Correlation of the relative fluctuations in squared density, rr , and temperature, T, and their lcr 
error bars. The solid lines give the adiabatic exponents of y = 5/3 (ideal gas) and y = 4/3. The third line 
shows the anticorrelation of density and pressure expected for contact discontinuities. The positive correlation 
not far from an ideal gas exponent indicates that adiabatic pressure fluctuations, as expected for the ram 
pressure effects of turbulent motions dominate also the density and temperature fluctuation spectrum. Right : 
Resulting shot noise subtracted, projected power spectrum of the pressure fluctuations observed in the central 
region of the Coma galaxy cluster (with 1<t errors). Also shown are power spectra for different power law 
shapes projected in an analoguous way as the observed spectrum. The model power spectra are labled with 
the exponents of the non-projected 3-dimensional power spectra. A Kolmogorov-Obuchov s pectrum with an 
expon ent of n = -7/3 is not far from the results, which lay in between n = 5/3 and n = 7/3 iSchuecker et al. 
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Fig. 20 Ratio of the 2p-3.s lines of Fexvii at 17.05&17.10 A to the 2p~3d line of the same ion at 15.01 A. 
The ratio has a significant peak at the center of the galaxy, NGC 4636. Since the lines originate from the 
same ion, the ratio cannot be due to a spatial abundance gradient. Given the narrow temperature range and 
the weak dependence of the line ratio on temperature it cannot be due to a temperature gradient either. This 
trend is the best evidence so far for resonant scattering of the line at 15.01 A in a hot galactic halo. 



the transition probabilities are large, photons get absorbed, but since the time between the 
absorption and emission is extremely short they get very quickly reemitted into a different 
direction. Because of the very short time between the abso rption and the emissio n of the 
photon the process can be effectively regarded as scattering. iGilfanov et alj dl987t) pointed 
out that because the optical depth of the resonance line depends on the characteristic velocity 
of small-scale gas motions, measurements of this optical depth give important information 
about the turbulent velocities in the hot plasma. 
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The first constrain ts on turbulent velocities using resonant scattring were obtained by 
IChurazov et alJ d2004) . who used XMM-Newton EPIC data to compare the relative fluxes 
of the He-like Fe K a and Fe lines in the core and in an annulus around the core of the 
Perseus cluster. Since the Fe K a at 6.7 keV has a much larger optical depth than the line 
a difference in their ratios in the two spatial regions with diff erent column densities would 
be an evidence for resonance scattering in the core of Perseus. Ichurazov et al.lll2004) found 
no evidence for resonance scattering in Perseus, indicating that differential gas motions on 
scales smaller than ~ 100 kpc in the core of the cluster must have a range o f velocities of at 
least half of the sound speed. Independently, Gastaldello & MolendJ ( 12004) reached similar 
conclusions. 

The first unambiguous evidence for resonant scattering was found using high-resolution 
spectra of the hot halo aroun d the giant elliptical galaxy NGC 4636 obtained by XMM- 
Newton RGS dXu et al J 12002). The plasma with temperatures below 0.9 keV observed in 
elliptical galaxies, groups of galaxies, and in the cool cores of some clusters emits three 
strong Fexvn lines. While the line at 15.01 A has a very strong oscillator strength and is 
expected to be optically thick, the blend of lines at 17.05 and 17.1 A have negligible optical 
depths. The radial profile of the ratios of these two lines (Fe xvn 17. 1A / Fe xvn 15.0 A) de- 
rived using RGS shows a clear gradient with the peak in the center of the galaxy NGC 4636, 
proving that m any of the 15.01 A photons get scattered before exiting its core (see Fig. [20] 
IXu et al.ll2002h . 

Following up the work by IXu et al.1 ( l2002h . IWerner et all J2009h analyzed the XMM- 
Newton RGS data of five nearby bright elliptical galaxies and found that the Fe xvn lines in 
the cores of four galaxies show evidence for resonance scattering in the innermost region. 
The data for NGC 4636 in pa rticular allowed the effects of resonant scattering to be studied 
in detail. [Werner et af (2009) used deprojected density and temperature profiles obtained by 
Chandra to model the radial intensity profiles of the strongest resonance lines, accounting 
for the effects of resonant scattering, for different values of of the characteristic turbulent 
velocity. Comparing the model to the data they found that the isotropic turbulent velocities 
on spatial scales smaller than «1 kpc are less than 100 km s" 1 and the turbulent pressure 
support in the galaxy core is smaller than 5% of the thermal pressure at the 90% confidence 
level, and less than 20% at 95% confidence. Note that the spatial scales of turbulence probed 
in the cores of elliptic al galaxie s by RGS are much sm aller than those probed in Perseus 
dChurazov et al.ll2004l) or Coma dSchuecker et alj|2004h . 

The Fe xvn lines are the most sensitive probes of turbulence using resonant scattering, 
but because they are present only at temperatures less than ~1 keV they cannot be used to 
probe turbulence in the much hotter clusters of galaxies. High-resolution spectra obtained 
by X-ray calorimeters on the future satellites like Astro-H and 1XO will allow us to probe 
turbulence also in higher mass systems and at larger radii (see Sect. [7}. 



4 AGN Heated Cluster Cool Cores 

The temperature profiles shown in the right hand panel of Fig.[l0]reveal two types of clus- 
ters: those with temperature profiles falling towards the center are called cool core (CC) 
clusters; and clusters with no central drop of temperature are called non-cool core, NCC 
clusters. While the former also have high densities in their cores which implies short cooling 
times for the central ICM, typically one or two orders of magnitude smaller than the Hub- 
ble time, the NCC clusters have usually central cooling times exceeding the Hubble time. 
This observation implies that the ICM plasma should cool and condense in the CC clusters 
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Fig. 21 Sketch of the temperature evolution of a clas- 
sical cooling flow model, with a reservoir of only 
slowly cooling hot gas at low density and a steady 
cooling of the fraction of gas with a mass cooling 
rate, M, down to lower temperatures. The spectrum 
of this cooling flow is then determined by the inte- 
gral of the spectra of the transient temperature phases, 
where each temperature phase contribution is fixed by 
the amount of heat that has to be radiated in the given 
temperature interval to cool an amount of ICM with 
mass M (see equations in the text). 



in the absence of any heat source, which could balance the cooling. Theoretical considera- 
tions focussed on the consequences of cooling withou t heating lead to the so-called "cooling 
flow" scenario (Fabi an & Nulsedll977uFabianlll994l) . which predicts mass deposition in the 
centers of some clusters at a rate of up to hundreds to thousands of solar masses a year. In 
general mo re than half of the clusters in the nearby Universe should have cooling flows in 
this model dPeres et al.lll998h . where the exact fractional number depends to some degree 
on the definition of a cooling flow. 

Interestingly, in particular in the context of this review, the cooling flow model leads 
to the prediction of very specific and testable spectral signatures of a steady state cooling 
flow. Looking at the plasma distribution in this scenario in temperature space, we have a 
reservoir of hot plasma and a steady cooling of some plasma with a mass deposition rate 
M, as sketched in Fig. 21. This simplified picture of a reservoir of hot gas and a constant 
fraction of gas cooling is motivated by the fact that cooling is accelerating with decreasing 
temperature because the plasma is getting denser while it is cooling in pressure equilibrium 
and the cooling radiation power is proportional to the gas density squared. For the part of 
the plasma that is steadily cooling to lower temperature, each temperature interval has to 
contribute to the cooling radiation with a power equal to the enthalpy change of the plasma 
that is cooling across this temperature interval: 

L x dT = —Enthalpy dT = 5 ^ dT m. (1 1) 

dt 2/jnip 

To calculate the expected spectrum we have to weight this power by the ratio of the emission 
at a certain frequency to the bolometric emission: 

emissivity(vWv 5k h dT . 

LJv) dvdT = f^— x — M. (12) 

emissivity io/ 2fim p 

To obtain the total spectrum radiated by all steady cooling temperature phases we have to 
integrate the above equation over dT: 

5k b dT . f r '»» A v (T')dv 

L x (v)dv = M - — -f-dT , (13) 

2/im p JT cmft A bo i(T') 
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Fig. 22 XMM-Newton RGS spectr um of the cen tral region of the prominent cool core cluster, A1835 com- 
pared to model predictions I Peterson et al. 2001). The red line shows a spectral model for an isothermal 
temperature of 8.2 keV. The blue line shows the expectation for the classical cooling flow and an ambient 
temperature of 8.2 keV, the green line shows a cooling flow model where temperatures below 2.7 keV have 
been truncated. Clearly, some of the cooling flow predicted lines are missing from the observed spectrum. 

where A V (T) and Aboi(T) are the emissivity of the plasma (radiation power per unit emission 
measure) for radiation at frequency, v, and for bolometric radiation, respectively. 

With the advent of high quality spectra provided by XMM-Newton and Chandra, the 
cooling flow paradigm could be tested by spectroscopy. First evidence that the cluster cool 
cores are not consistent with the classical cooling flow model came from spectra analyzed 
with the Reflection Grating Spectrometers on XMM-Newton. Some of the cool core clusters 
have such highly peaked surface brightness profiles, that they look almost like blurred point 
sources to the RGS (which operates analogously to slitless spectroscopy in optical astron- 
omy). Therefore the RGS can obtain spectra with much higher resolution than provided by 
the imaging CCD devices of the XMM-Newton EPIC cameras for these very peaked cool 
core regions. Fig.l22lshows the first published R GS spectrum of one o f the very strong cool 
core clusters, A1835, at a redshift of z = 0.2523 jPetersone t al. 2001). In the Figure the ob- 
served XMM-Newton spectrum (dark red) is compared to the prediction for the cooling flow 
model with a mass deposition rate inferred from the total radiative emission in the cooling 
flow region determined from the X-ray images (blue curve). There are clearly some expected 
lines, notably those of Fexvn, missing in the observations. A more detailed analysis of the 
spectrum shows that for an ambient ICM temperature of 8.2 keV, the central regions shows 
temperature phases down to about 2.7 keV, but the expected intermediate temperature phases 
below 2.7 keV are absent in the spectrum. This was generally taken as the first evidence that 
there is something wrong with the classical cooling flow model. Since at the same time high 
angular resolution Chandra images showed signs of strong interaction of the central AGN 
with the ambient I CM in many cool core regions (e.g. lDavid et al ]|200lLlNulsen etalj2002l ; 
iFabian et alj|2005l) . the most probable solution to the absence of strong cooling was soon 
believed to be the heating of cool core regions by the central AGN. 

The spectral signature of the absence of cooling was not only obtained from XMM- 
Newton RGS spectra. For example for the cool core region in the X-ray halo of M87, the 
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Fig. 23 Left: Spectrum of hot plasma at various temperatures around the Fe L-shell line blend near 1 keV 
as observed with the response function of the XMM-Newton EPIC pn instrument. This feature lends itself as 
a nice thermometer for the temperature range 0.2 to 2 keV. Middle: Predicted X-ray spectrum around the Fe 
L-shell line blend for the cooling flow in the X-ray halo of M87 with a mass deposition rate of M ~ 1M yr~' 
compared to the actual spectrum in the radial region r = 1 -2 arcmin as observed with the XMM-Newton EPIC 
pn instrument. The observations are clearly inconsistent with the predictions from the classical cooling flow 
model. Right: Reproduction of the observed spectrum by assuming that the plasma in this re gion of the M87 
halo occupies the narrow temperature interval of about 1.44 to 2 keV (for more details see Bohringer et al. 
[2002b . 



central dominant galaxy in the northern part of the Virgo cluster, clear evidence for an in- 
consistency of the spectral signatures with a classical cooling flow model could also be 
provided by results with the EPIC pn and MOS cameras of XMM-Newton. The most promi- 
nent feature that can be used for these diagnostics is the blend of Fe lines in the spectrum 
from transitions into the L-shell at photon energies around 1 keV. Fig. [23] shows how the 
spectral feature of this blend of Fe lines changes with temperature. The shift is caused by 
the changing degree of ionization of the Fe ions that contribute to the line blend. When the 
temperature is lowered more and more electrons in the Fe ions shield the charge of the nu- 
cleus decreasing the binding energy of the L-shell electrons. As illustrated in Fig. 123 1 this 
spect ral feature can be use d as a sensitive thermometer in the temperature range of 0.2-2 
keV (Boh ringer et al ]2002l) . The middle panel in Fig.l23lshows the spectrum implied for the 
broad temperature distribution expected for the cooling flow in M87 (in the radial range 1-2 
arcmin) with a mass deposition rate of M ~ 1M yr" 1 . The predicted spectrum is clearly 
different from the observed one. A way to modify the predicted spectrum in the direction 
of the observed one, is to assume that the central region of M87 is highly absorbed by cold 
or warm material to reduce the low energy part of the broad Fe line feature. The actual ab- 
sorption can be determined directly, however, by using the power law spectrum of the M87 
AGN with the result that there is no significant excess absorption above the known value 
of the galactic absorption column density of about 1.8 x 10 20 cm -2 . A very good fit to the 
observed spectrum is then obtained by assuming that the bulk of the plasma in the central 
region of the M87 halo only occupies a narrow temperature interval of about 1 .44 to 2 keV 
(Boh ringer et al.ll2002h as shown in the right panel of Fig. [23] 

It has to be noted that the interpretation of spectra from cooling flows had an interesting 
history alr eady before t he lau nch of Chandra and XMM-Newton. For example llkebe et all 
l i 19971) and lMakishimal i ll 999b already argued on the basis of ASCA spectroscopy with simi- 
lar arguments as those given above for M87, that the prediction of cooling flow models is not 
met by the observations. The counterargument to keep the cooling fl ow model alive at the 
time was to postulate internal absorption (see e.g. lAllen et al.ll200ll) . The spectral resolu- 
tion and sensitivity of ASCA was not quite good enough, however, to reject the cooling flow 
model in the very compelling way as now done with the XMM-Newton results. Even earlier, 
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Fig. 24 Left: Differential emission measure distribution as a function of temperature of the plasma in the cool 
core region of the cluster Abell 262 as deduced from a spectral analysis of XMM-Newton RGS observations. 
The results are derived for the innermost four shells: - 0.5 armin (triangles), 0.5 - 1 armin (squares), 1-2 
arcmin (stars), 2-3 arcmin (circl es). For comparison the dashed line shows an approximation to the slope of 
the isobaric cooling flow model (Kaastra et al. 2004). Middle: Differential emission measure distribution for 
the isobaric cooling flow model, for 0.5 solar abundance I Kaastra et al. 2004) Right: Differential emission 
measure distribution of the plasma in the inner 2 arcmin region in the core of the Centaurus cluster of galaxies. 
The results from a spectral analysis of XMM-Newton RGS and CHADRA are shown for comparison. The 
solid line shows the expected emission measure distribution f or a 10 M Q yr -1 isob aric cooling flow at solar 
metallicity, cooling from a temperatue of 4.5 keV to 0.08 keV ( Sanders et al. 20081). 



spectroscopic analysis of the X-ray halo of M87 with the Crystal Spectrometer on board 
of the Einstein satellite showed th e emission lines for the cool phases expected in cooling 
flows JCanizares et al ][l97i,[l982), that we are now missing in modern observations. This 
evidence for steady state cooling in cooling flows gave actually one of the strongest supports 
to the cooling flow scenario at the time. It is now clear, however, that the results obtained in 
these early days were most probably features of spectral noise. 

Deep observations of cool core regions, in particular with the XMM-Newton RGS instru- 
ment, have now also been used to characterize the multi-temperature structure in cool cores 
in more detail in terms of the emi ssion measure (or luminosity) distribut ion of the plasma 
as a function of temperature (e.g. lPeterson et a l. 2003; Kaastra et all2004T) . This works well 
only in colder systems. Fig.[24]shows in the left panel the differential emission measure dis- 
tribution in the cool core region of the cluster Abell 262 in the innermost four shells out to 3 
arcmin. These observations can be compared to the expected emission measu re distribution 
for th e model of an isobaric cooling flow, shown in the middle panel of Fig.|24]( Kaastra et al. 
120041) . While the temperature distribution is not isothermal, we clearly note that the emis- 
sion measure distribution falls off much more steeply than the model expectations towards 
lower temperatures, highlighting again the fact that less plasma is observed at lower temper- 
atures than required for the cooling flow model. A similar result is shown - in the righ t panel 
of the Figure - for the case of the Centaurus cluster taken from I Sanders et aD §008) based 
on a very deep XMM-Newton observations of the cluster with the RGS instruments as well 
as deep Chandra data. Here we can follow the local temperature distribution of the plasma 
over a full decade in temperatures from about 0.4 to 4 keV. But again the contribution of the 
very low temperature phases is much less than what is expected for unimpeded cooling. 

As the change from the cooling flow paradigm to an AGN feedback scenario seems to be 
now widely accepted, the interest has shifted to the question: how is the ICM actually heated 
by AGN interaction? One of the most important earlier arguments in favour of cooling flows 
in the absence of feedback was that any possible heating mechanism has to be very well fine 
tuned, to exactly provide the balance to cooling. If too much heat is produced, it disperses the 
observed dense gas cores in cooling flows. Secondly, we observe increasing entropy profiles 
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Fig. 25 False colour image of the central region of the Perseus cluster produced from images in three energy 
bands 0.3-1.2 keV, 1.2-2 keV, and 2-7 keV. An image smoothed on a scale of 10 arcsec (with 80% normal- 
ization) has been subtracted from the image to highlight regions of strong density contrast. The image shows 
a series of nearly concentric "rippl es" which are interpreted as sound waves or weak shock waves set off by 
the activity of the central AGN iFabian et all2006l) . 



in the ICM of clusters (Fig. \T3[ . It therefore has to be explained how central heating can 
work without inverting the entropy profiles leading to the dissipa tion of heat b y convection. 
Both requirements need to be met and seem almost improbable (Fabian| [l994h . This makes 
it even more interesting now to understand how nature achieves this fine tuning. 

One of the first, best studied cases is that of the Perseus cluster, which has now been 
observed wi t h the Chandra observatory for more than 1 Million seconds as proposed by 
IFabian et alJd2005h providing the most detailed picture of a central cool core re gion. Fi g. [25] 
shows a multi-colour image of the central region of the Perseus cluster from IFabian et al.l 
J2005I) . The "false colour" image was produced from images in three energy bands 0.3- 
1.2 keV, 1.2-2 keV, and 2-7 keV. An image smoothed on a scale of 10 arcsec (with 80% 
normalization) has been subtracted from the image to highlight regions of strong density 
contrast. The central, bright part of the image shows two features of low surface brightness, 
interpreted as two cavities of the X-ray emitting plasma filled with the relativistic plasma 
from the radio jets of the central AGN. This region is surrounded by a series of nearly 
concentric "ripples" which are more clearly brought out by the unsharp masking processing 
of the image. 

A more detailed analysis of the density and temperature variations across the ripples, 
shown in Fig. [26] implies typical pressure variations associated with the ripples with ampli- 
tudes of about ±5-10%, which are i nterpreted as sound waves or very weak shock waves in 
the innermost region. As detailed in IFabian et al l (T2005) the dissipation of the sound waves 
- if the viscosity is sufficiently high - can balance the radiative cooling in the cool core by 
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Fig. 26 Temperature, density, pres sure, and pre s sure r esidual profiles in the core region ICM of the Perseus 
cluster in a sector to the North-East Fabian et al. ( 2006). The red line shows the ripples from an unsharp mask 
image. The dashed line marks the position of a shock front. 



heating with the realistic assumption that the sound waves cross the cool core region in about 
1/50 of the cooling time. 

iFabianelail 12006) have also carried out multi-temperature plasma modeling in various 
regions of the cluster centre and produced mass distribution maps of plasma at temperatures 
around 0.5, 1, 2 and 4 ke V. These maps are sho wn in comparison to a map of the optical 
Ha emission (observed bv lConselice et al.ll200ll) in Fig. [27] The distribution of the lowest 
temperature phase around 0.5 keV has a striking similarity to the Ho- map. There is relatively 
little plasma at 1 keV, mostly in the very center of the cool core. The 2 keV phase coincides 
mostly with the denser regions of the ICM. This detailed picture of the central region of the 
Perseus cluster given in a series of papers by Fabian et al. and Sanders et al. well illustrate 
the frontier of the application of X-ray spectroscopy and thermodynamic analysis of the 
ICM in the bright central regions of clusters. 

To show that the physics of the cool core regions seems to be quite diverse and that we 
can obviously not easily generalize the behaviour of the Perseus cluster to all cool core clus- 
ters we show two further examples of systems studied in detail, that of M87 and of Hydra 
A. M 87 was studied in very m uch detail by a combination of a 500 ks ex posure with Chan- 
dra IIForman et al .l2005l,l2007h and a 120 ks exposure with XMM-Newton dSimionescu et alj 
I2007L 2008b). The X-ray halo of M87 is the closest of all cool core regions in galaxy clus- 
ters and can therefore be studied at the highest spatial resolution. Also the AGN in the 
center of M87 is active and is interacting wit h the ambient ICM through relativistic jets (e.g. 
iBohringer et alll 19951 : IChurazov et a"i]|200ll) . The deep Chandra exposure shows a surface 
brightness discontinuity which can be almost followed over the entire circumfer ence of a 
ring around the nucleus with a radius of 14 kpc (Fig. |28lfForman et alj2 005. 2007). Another 
arc-shaped discontinuity is seen over a s maller sector region a t the radius of 17 kpc. T hese 
features are identified as shock waves in iForman et af] (2005) and lForman et al. 1 120071) . A 
detailed spectroscopic analysis of the M87 X-ray halo with XMM-Newton can corroborate 
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Fig. 27 Ha image in the halo of NGC 1275 in the center of the Perseus cluster (from Conselice et al. 12001), 
upper left) shown on the same scale for comparison with maps of the mass distribution of different tem- 
peratu re components, at 0.5, 1, 2, and 4 keV, of the ICM in the center of the Perseus cluster I Fabian et al. 



this interpretation and give enhanced physical constraints. While the density discontinuity 
seen with Chandra implies a shock Mach number of about 1.2, the temperature jump seen 
in the projected spectra is about 5% at about 2 keV and implie s a shock Mach number of 
> 1.05 without deprojection correction dSimionescu et al] |2007h . Thus the temperature en- 
hancement identifies the surface brightness discontinuity definitely as a weak shock front. 
The shock front has also been visualized in a pressure map of the M 87 halo constructed 
from the XMM-Newton data shown in Fig. 1281 ( Si mionescu et al]|2007h . Similar to the hy- 
drodynamic maps described above this map has been derived from spectral modeling in 
many spatial pixel, which in t his case have been con structed by means of a Voronoi tessel- 
lation technique developed by Die hl & Statleri < t2006h and originally by Cappell ari & Copinl 
d2003h such that each pixel has a spectral signal to noise of at least 100. To enhance the 
visibility of discontinuities in the presence of the strong pressure gradient, a smooth, spheri- 
cally symmetric model has been subtracted from the pressure distribution shown in Fig. [28] 
The map clearly shows a pressure enhancement at the shock radius of 3' (14 kpc) which im- 
plies without projection correction a Mach number larger than 1.08. There are some further 
interesting features seen in the pressure map. Outside the shock front we note two regions 
of lower pressure in the NE and in the SW. These are the regions where radio lobes filled 
with relativistic plasma emerging from the central AGN are observed. This observation thus 
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Fig. 28 Left: Chandra image of the X-ray halo of M87 in the 0.5-2.5 keV band iForman et al.l2005l) . Several 
features are marked on this image: the surface brightness enhancements coinciding with the radio lobes (SW 
arm and E arm) and several sharp surface brightness discontinuities (14 kpc ring, 17 kpc arc, and 37 kpc 
arc) which are identified as shock fronts. Right: Map of the pressure deviations from a smooth spherica l 
symmetric model derived from a deep XMM-Newton observation of the M87 halo iSimionescu et al. 2007). 
The 14 kpc ring is now clearly revealed as a pressure discontinuity supporting the shock interpretation. The 
regions of the radio lobes (SW arm and E arm) are now showing up as zones of pressure deficit, which is 
probably compensated for by the unseen pressure of the relativistic radio plasma. 



implies that (at least part of) the missing pressure in these regions comes from the pressure 
of the relativistic plasma in the lobes, which does not contribute to the X-ray emission anal- 
ysed here. Further insight into the physics of this interesting cool core region in M87 will be 
described in connection with the analysis of the element abundance distribution discussed 
in the next section. 

In a combined Chandra and XMM-Newton study of the central region of the galaxy 
cluster associ ated with the radio sour ce Hydra A, a similar fron t diagnostics has also been 
performed by Nuls en et alj d200 §t) and lSimionescu et"ai] ( 120083) . Here the shock feature has 
a radial extent of 200 - 300 kpc and hints at an event which was much more powerful by 
about two orders of magnitude than that in M87. Fig.|29]shows again the surface brightness 
discontinuity cle arly detected in the high angular resolution image obtained with the Chan- 
dra observatory dNulsen et al.l2 005). A deep XMM-Newton observation then provides higher 
photon statistics that allows us to study this discontinuity spectroscopically. Firstly, we can 
again visualize th e shock front as a pressu re jump in the XMM-Newton based pressure map 
shown in Fig. [29l(S imionescu et al.l2008ah . A more detailed analysis of the temperature and 
density distribution around the shock front in several sectors, shown in Fig. [30] provides 
more detailed information about the shock strength. Modeling the discontinuity in three di- 
mensions and fitting the proj ected model to the data i mplies a Mach number of the best 
fitting model of about 1.2-1.3 dSimionescu et al . 2008a). Simple one-dimensional modeling 
of the evolution of the shock driven by the energy input of the radio lobes as well as detailed 
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Fig. 29 Left: Chandra image of Hydra A with the radio lobes superposed (Nulsenetal. 2005). An oval 
shaped surface brightness discontinuity at a radius of about 200 - 300 kpc can be noted that is identified with 
a shock wave caused by an energetic radio jet outburst. Right: Map of the pressure deviati ons from a smooth 
azimu fhally symmetric model derived from a deep XMM-Newton observation of Hydra A I Simionescu et al. 
2008a). The location of the X-ray surface brightness discontinuity observed in the Chandra image is approx- 
imated by the elliptical black line. A clear pressure jump can be seen almost all along this line indicating a 
shock front with a Mach number of about 1 .3. 




Fig. 30 Left: Model of the effect of a shock front with different Mach numbers predicting the o bserved 
surface brightness distribution in the XMM-Newton image of Hydra A I Simionescu et al. 2008a) Right: 
Temperature profile around the shock front compared to a smooth model (blue lines). The observed tem- 
perature and density (su rface brightness) jumps are consistent with a shock wave with Mach number ~ 1.3 
(Simionescu et al. 2008a). 



36 




Energy (keV) 



20 30 
Wavelength (A) 



Fig. 3 1 Left panel: The li ne spectrum of the cluster 2A 0335+096, as observed with XMM-Newton EPIC 
(from Werner et al. 2006b). Right panel: Line spectrum of M 87, as observed with XMM-Newton RGS (from 
Werner et al. 2006a). This representation of the spectrum has been obtained by recalculating the spectrum 
with the best fit values without lines, and using the result to subtract it from the observed spectrum as indicated 
on the Y-axis. While the CCD spectra are plotted as a function of the observed energy, the grating spectra are 
shown as a function of the observed wavelength. 



three-dimensional hydrodynamical simulations provide good estimates for the age and the 
total energy of the eve nt causing the shock wave, with values for th e age of ~ 2 x 10 7 yr and 
the energy of 10 61 erg dNulsen et al.ll2005l ; ISimionescu et alj|2008cj) . 

These few examples show the complex and probably very diverse physics prevailing in 
cluster cool core regions. We hope that further detailed studies of a larger sample of cool core 
regions will help to generalize the scenario of cool core heating. A major break through is 
also expected here from future high resolution spectroscopy, which will allow us to observe 
turbulent and bulk motions in the ICM as described in section 7. 



5 What Chemical Abundance Measurements Tell Us 

Clusters of galaxies are unique laboratories for the study of the nucleosynthesis and chemical 
enrichment of the Universe. Their deep gravitational potential wells keep all the metals 
produced by the stellar populations of the member galaxies within the cluster. The dominant 
fraction of these metals reside within the hot ICM. The chemical abundances measured in 
the intra-cluster plasma thus provide us with a "fossil" record of the integral yield of all 
the different stars (releasing metals in supernova explosions and winds) that have left their 
specific abundance patterns in the gas prior and during cluster evolution. 

X-ray spectroscopy provides an accurate measure of metal abundances in the ICM. 
These abundances put constrains on nucleosynthesis and on the star formation history of 
the clusters, and as long as the stellar populations where the cluster metals were synthesized 
are representative, these measurements allow us to put constraints on the chemical evolution 
of the Universe as a whole. 

The CCD type detectors on Chandra, XMM-Newton, and Suzaku allow us to detect the 
emission lines of O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni (see the left panel of Fig. 13 lb. The 2p- 
ls Ne lines at 1.02 keV are in the middle of the Fe-L complex (lying between 0.8 to 1.4 keV). 
The energy resolution of the CCD type detectors is not sufficient to resolve the individual 
lines in the Fe-L complex, which makes the Ne abundance determination very sensitive to 
uncertainties in the temperature modeling. However, the abundance of Ne in the cores of 
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cooling core clusters, observed with XMM-Newton with sufficiently deep exposures, can be 
well determined using the high spectral resolution of the Reflection Grating Spectrometers 
(RGS). The RGS has a higher spectral resolution than the CCDs. It resolves the Fe-L com- 
plex into individual lines and allows for accurate measurements of the equivalent widths of 
O, Ne, Mg, and Fe. In the case of deep observations of nearby bright, relatively cool clusters 
and elliptical galaxies even the spectral lines of C and N can be detected (see the right panel 
of Fig. [3Tb. The high resolution transmission grating spectrometers on Chandra are because 
of their small effective area unfortunately not well suited for cluster spectroscopy. 

The equivalent widths of the observed lines are under the assumption of collisional equi- 
librium directly converted into abundances of the corresponding elements. Because the ICM 
is a relatively uncomplicated physical environment (ionization non-equilibrium or optical 
depth effects are minimal) these abundances are relatively robust. 

Abun dances are often shown in the literature with respect to the outdated Solar abun- 
dances by Anders & Grevesse ( 1989) or with respect to newer sets of Solar and proto-solar 
abundances by Lodders (2003). The more recent Solar abunda nce determinations of O, Ne , 
andFeby Lodders (2003) are ~30% lower than those given bvl Anders & Grevessd ( fl989h . 



The Solar abundances of O and Ne reported by iGrevesse & Sauval 
the abundance of Fe is slightly lower than those reported bv lLodders 



1998) are higher and 



2003). 



5.1 Early pioneering work on chemical abundances 

The first evidence that the ICM is strongly polluted by metals ejected from stars in the 
cluster galaxies came with the discovery of the Fe-K line emission in the spectrum of the 
Perseus, Coma, and Vir go clusters by the Ariel V and OSO-8 satellites iMitchell et al.l 19761 ; 
ISerlemitsos et alj|l977t) . Spectral analysis of samples of clusters observed with OSO-8 and 
HEAO-1 A2 showed that the ICM is enric hed by Fe to 0.3-0.5 of the Solar abundance value 
jMushotzkv et al.lll978l ; [Mushotzkvlfl98i) . 

Until the launch of ASCA in 1993, only the Fe abundance was accurately measured in 
many clusters. ASCA for the first time detected emission features from many different el- 
ements (O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni) in the ICM, opening thus a new chapter 
in the chemical enrichment studies. While core collapse supernovae (SNcc) produce large 
amounts of O, Ne, and Mg, Type la supernovae (SN la) produce large quantities of Fe, 
Ni, and Si-group elements (Si, S, Ar, and Ca), but only very little O, Ne, and Mg. Early 
ASCA measurements of O, Ne, Si, S, and F e abundances in four c lusters suggested that 
the ICM enrichment is dom inated by SNcc ( Mushotz kv et al.ll 1996b . However, few years 
later analyzing ASCA data, iFinoguenov et alj ( 12000b showed that while the SNcc prod- 
ucts are uniformly distributed in the ICM, the chemical enri chment of the cluster c ores 
is dominated by SN la. A n interesting result based on ASCA by Fuka zawa et alj d 1998b and 
iBaumgartner et all fe005b showed that the Si abundance and the Si/Fe ratio increase from 
the poorer to the richer clusters, suggesting that the relati ve contribution of SNcc increases 
toward the more massive clusters. In a pioneering work iDupke & Whitel ( 12000) used the 
abundances measured in three clusters of galaxies to put constraints on SN la models. They 
used the Ni/Fe ratio to discriminate between different SN la explosion scenarios available 
in the literature. Their work preferred the W7 "deflagration" models over the "delayed det- 
onation" models. 

The ASCA data of th e Centaurus cluster we r e the first to show a str ong Fe abundance 
peak at the cooling core (Fuka zawa etafl 1 1994b . lAllen & Fabian] J 1998b used ASCA spec- 
tra and ROSAT images to show a segregation between the metallicities of clusters with and 
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Fig. 32 The average radial distribution of the Fe abundance for a sample of cooling core (filled circles) and 
non-cooling core (empty circles) clusters observed with BeppoSAX (from De Grandi et al. 2004). 



without cooling cores and suggested that it is caused by the presence of metallicity gradients 
coupled with the sharply peaked surface brightness profiles of the cooling core clusters. Un- 
fortunately, the large point-spread function of ASCA did not allow to investigate the spatial 
abundance distributio ns in detail. Using the BeppoS AX satellite, which had a better spatial 
resolution than ASCA, |Pe Grandi & Molendi (2001) clearly showed that while non-cooling 
core clusters have flat Fe abundance profiles, the metallicity distribution of cooling core clus- 
ters has a gradient with a central peak (see Fig. [32l. Abundanc e gradients were also found 
in 9 out of 10 galaxies and groups of galaxies bv lBuotd ( I2000a[) . who analyzed ROSAT data 
of these ~ 1 keV systems. 

ASCA data were furthermore used for a pioneering search fo r evolution in cluster metal- 
licities, revealing a lac k of evolution out to redshift z ~ 0.4 dMushotzkv & Loewensteir] 



1997 : lRizza et al.l l998) and finding no evidence for a decrease at higher redshifts jDonahue et al.l 
19991) . 



5.2 Chemical enrichment in cooling core clusters 

The emerging picture of the chemical enrichment in clusters of galaxies at the time of the 
launch of XMM-Newton and Chandra was that of an early enrichment of the ICM by SNcc, 
the products of which are today well mixed and homogeneously distributed, and a subse- 
quent more centrally peaked enrichment by SN la which continue to explode in the cD 
galaxy for a long time after the cluster is formed. According to this picture the observed 
peaked distribution of the Fe abundance is largely due to the SN la in the cD galaxy. 

The first abundance studies performed with XMM-Newton seemed to confirm this sce- 
nario. The radial abundance profiles of the cooling core cluster Abell 496 showed a flat dis- 
tribution of O, Ne, Mg, which are pre dominantly produce d by SNcc, but a centrally peaked 
distribution of Si, S, Ar, Fe, and Ni ( iTamura et al.ll200l[) . However, while most of the Fe 
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and Ni is produced by SN la, Si, S, Ar, and Ca (so called a elements) are produced by both 
supernova types in similar proportions. In many systems, Si has a relatively well determined 
abundance value and in a cluster with a flat O profile one would expect to see a shallower Si 
abundance gradient than that of Fe. 

However, this has never been observed. In the most nearby cluster cooling core, M87, the 
Si abundance has a similar gradient as that of Fe, b ut at the same time it has a low O/Fe ratio 
in the core that increases toward the outer regions (Bohringer et al. 200ll: lFinoguenov et al.l 
|2002|) . A similar abundance pat tern was observed in the cluster sa mple analyzed bv | Tamura et al.l 
J2004I) , in the Per seus cluster dSanders et alj|2004h. in Abell 85 tourret et al.ll2005[) . in the 
Centaurus cluster dlvlatsushita et alJ2007al:ISanders & Fabianl200d) , in the group NG C 5044 
feuote et al.ll2003l) , in Abell 1060 jsato et alj|2007h. and in AW M7 dSato et al.ll2008l) . 

To solve this discrepancy in M 87, lFinoguenov et "aTl d2002h proposed that there are two 
types of SN la with different Si/Fe yields. Those with longer delay times and higher Si 
yields are today dominating the enrichment of elliptical galaxies, such as the cD galaxies in 
the centres of cooling core clusters, and those with shorter delay times and lower Si yields 
dominate the enrichment in the ICM at larger radii, because the bulk of them exploded in 
the time of the strongest star formation period (for more details on SN la models see the 
Subsect.[53]l. 

Recently, ISimionescu et al ] d2009h determined radial abundance profiles for O, Si, S, 
and Fe in the Hydra A cluster. They found that all elements have a centrally peaked distribu- 
tion. In order to investigate the radial abundance distributions in clusters in more detail, they 
combined their profiles with those determined for 5 other clusters with high quality X-ray 
data (M87, Centaurus, Fornax, 2A 0335+096, Sersic 159-03, and Abell 1060). In the com- 
bined data set they find that the O distribution also peaks in the cluster core and decreases 
with radius as dO/d(\og r/r2oo)=~0.48±0.07. But the Fe abundance still decreases faster: 
dFe/d(\og r/r 2 oo)=-0.72±0.04, and therefore the fraction of O with respect to Fe increases 
with radius as d(0/Fe)/d(\og r/r2oo)=0.25±0.09. The Si abundance profile, however, remains 
exactly the same as that of Fe, it decreases with radiu s as: dS i/d(log r/r2oo)=-0.72±0.04. 

Surprisingly, for the group NGC 507, 1 Sato et all d2009) found steep abundance gradi- 
ents peaking on the cluster core not only for Si, S, and Fe but also for Mg. They also found 
a mildly centrally peaked abundance distribution for O. The Suzaku satellite has a better 
sensitivity in the O line energy and it allows for better measurements of the Mg lines in 
the low surface brightness cluster outskirts than XMM-Newton, which has strong instrumen- 
tal lines and effective area calibration problems at the Mg energy. As shown in Fig. [33] the 
Mg/ Fe ratio as measured with Suza ku is consistent w ith being flat out to 0.2ri8o in Forna x 
^Matsushita etai1l2007al) . HCG 62 ( iTokoi et al.|[2008h . and in NGC 507 JSato et al.|[2009h . 
Since both Mg and O are produced by SNcc their radial profiles should follow each other. 
However, even though the error-bars are large, the O/Fe ratio seems to increase with the 
radius in these three systems. The different abundance profiles for Mg and O might be due 
to systematic uncertainties in regions with low surface brightness where the confusion with 
the strong O vm Galactic foreground line becomes an issue. The most straightforward inter- 
pretation of these profiles is that the SN C c products are also centrally peaked, as indicated 
by Mg, and the measured peak in the O abundance profile is shallower because of the con- 
fusion with the Galactic O vm foreground lines which become important in the faint cluster 
outskirts. 

Based on the Fe abundance profiles, foe Grandi et alj J2004I) and lBohringer et all ( |2004|) 
concluded that the central abundance peak is due to the metals released by the stellar popu- 
lation of the cD gal axy. To produce the observed abundance peaks requires long (>5 Gyr) 
enrichment times (Bohr inger et alj |2004). In the absence of mixing the metallicity profiles 
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Fig. 33 Radial pro file of abundance ratios for NGC 507 iSato et alj|2009l) . HCG 62 iTokoi et all2008t) . and 
the Fornax cluster I Matsushita et al. 2007a). All the abundance ratios are consistent with having a flat radial 
profile, except O /Fe which shows an indica tion for an increase with the radius. The abundance ratios are in 
the solar units of Anders & Grevesse 1 1989). 



should follow the optical light profiles of the cD galaxies. But the observed metal mass 
profiles are much less peaked than the optical light, which implies that the metals get 
mixed and transported out to larger radii, most probably by the AGN/ICM interaction (e.g. 
iRebusco et al.ll2005ll200fjh . 

The role of the AGN in the metal transport was both theoretically and observationally 
best investigated in M 87. IChurazov et"aT]j200lh proposed a scenario where buoyant bubbles 
of radio emitting relativistic plasma, produced by the jets of the AGN, uplift the cool metal 
rich gas from the central parts of the cD galaxy. This scenario was confi rmed by observa- 
tions o f cool metal rich gas, distributed along the radio lobes (see Fig.[34j S imionescu et alj 
2008b). The Fe abundance of the uplifted gas (a mixture of ICM, stellar mass loss, and SN la 
products) is ~2.2 times solar and its mass is « 5 x 10 8 M Q . It takes »30-l 10 Myr to produce 
the metals in the cool gas. The relative abundances of the uplifted gas are consistent with 
those of the ambient ICM, indicating that the dominant fraction of metals in the cluster core 
was uplifted after the last major epoch of star formation and the relative enrichment rate by 
SN la with respect to that by stellar winds remains fairly constant with time. Uplift of metal 
rich cool gas from the cD galaxy by buoyantly rising bubbles of radio emitting plasma was 
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Fig. 34 Map o f the Fe a bundance in M 87 in Solar units indicated under the color bar. Contours of the 90 cm 
radio emission I Owen et al. 2000) are overplotted. Beyond the expected radial grad ient, one clearly sees th e 
enhanced Fe abundance in the radio arms, especially within the Eastern arm (from Simionescu et al. 2008b). 



also observed in the Hydra A cluster of galaxies dSimionescu et al The total mass of 

the uplifted Fe in Hydra A is 1.7 x 10 7 M Q . 

Can the abundance peaks be produced by only stellar winds and SN la in the cD galaxy? 
To explain the centrally peaked distribution of O, Si, and S in the Hydra A cluster would re- 
quire the amount of metals produced by stellar winds to be 3-8 times higher than predicted 
by mod els or the metals produc ed by SNcc in the proto-cluster phase not to be mixed com- 
pletely dSimionescu et al. 2009). Centrally peaked distrib ution of all elements c an also be 
produced by ram-pressure stripping of cluster galaxies (e.g. lDomainko et al.l2 006). But both 
simulations and observations indicate that galaxies are getting stripped already at large dis- 
tances from cluster cores and this mechanism probably does not contribute significantly to 
the observed abundance peaks on the relatively small spatial scales around the cD galaxies. 

The radial metallicity gradients in cooling flow clusters often display an inversion in 
the center (e.g. ISanders & Fabianll2002h . While an apparent drop in metal abundances can 
be the result of an oversimplified model in the spectral analysis of the very center of a 



cooling core (e.g.lBuote & FabiarJI 19981 ; iMolendi & GastaldelrolboOll ; IWerner et af]|2006bl 



iMatsushita et al. [2007al . for more discussion about spectral fitting biases see Subsect. r5.5.U 



in some cases this inversion seems robust. In particular the metallicity in the Perseus cluster 
drops very clearly at radii smaller than 40 kpc and the central dip does not disappear when 
multi-temperature models and extra power-law components are introduced in the spectral 
fitting, or when projection effects and possible effects of resonant scattering are taken into 
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Fig. 35 A compilation of the morphological type of the host galaxy versus the decline rate of SNe la. 
NGC 2841 has spectral features of an elliptical and NGC 0632 has emission lines characteristic for late- 
type galaxies. The vert ical markers indicate the average decline rates for early- and late-type galaxies. From 
Gallagher etai] <2005l) . 



account dSandersetai]|2004 ISanders & Fabianl 120061 120071) . The metallicity dips in the 
cluster cores may be related to the AGN/ICM feedback. 



5.3 Constraining supernova models using clusters 

As we mentioned before, clusters of galaxies retain all the metals produced by their stellar 
populations within their gravitational potential. T he dominant fraction of these metals re- 
sides within the hot ICM (Finogue nov et al] l2003). All the elements from oxygen up to the 
iron group are primarily produced by supernovae. Their abundances within the ICM thus 
provide us with the integral yield of all the supernova explosions in the cluster. 

Since the launch of the ASCA satellite there have been efforts to use the abundance pat- 
terns in the ICM to put constraints on the theoretical supernova mo dels. These efforts were 
mai nly focused on con straining the Type la explosion mechanisms. iDupke & Whitel bOOOl) 
and IDupke & ArnauJ feOOll) used the Ni/Fe ratio measured in clusters by ASCA to distin- 
guish between the SN la models assuming slow deflagration and delayed detonation. While 
the deflagration explosion mechanism of SN la (represented by the W7 and W70 models 



in the li terature , Iwamoto et al .11 19991) predicts high Ni/Fe ratio of 2.18 to 3.22 in the Solar 
units of iGrevesse & Sauvall 1 1998 ), the delayed-detonation explosion scenarios predict sig- 



nificantly lower Ni/Fe ratios of 0.9-1.4 Solar. They are represented by the WDD1, WDD2, 
WDD3, CDD1, and CDD2 models in the literature, where the last digit indicates the density 
at which the flame velocity becomes supersonic (deflagration-to-detonation transition den- 
sity) in units of 10 7 g cm -3 ; the "C" and "W" refer to two different central densities (1.37x 
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and 2.12 x 10 9 g cm 3 , respectively) in the model at the onset of the thermonuclear runaway 
lllwamoto et alJll999h . ASCA measured relatively high Ni/Fe ratios which clearly favored 
the W7 SN la model l lDupke & WhitdbOOfilDupke & ArnaucfeoOll see also Sect. 5.5.1). 

The results based on the abundances measured by XMM-Newton, Chandra, and Suzaku 
are not so clear and depending on the relative abundances in the given cluster sometimes 
deflagration models and sometimes delayed detonation models are favored. 

An excellent example of the complicated situation is the Virgo cluster with M87 in its 
cooling core. The abundance patterns in M 87 cannot be described by any combination of 
SNcc and deflagration SN la models. The low Ni/Fe ratio and very high Si/Fe ratios clearly 
favor delayed detonation models. The high SN la yield of Si-group elements in M87 imply 
that the Si burning in SN la is incomplete, as predicted by the delayed-detonation models 
with lower density of deflagration to denotation transition (WDD1 model). The abundance 
patterns in the outer region of M 87 are on the other hand characteristic of deflagration su- 
pernova models or delayed detonat i on m odels with the highest deflagration to detonation 
density (WDD3). iFinoguenov et al.l J2002h conclude that these abundance patterns confirm 
the diversity of SN la. Today, this conclusion seems to be supported by both optical super- 
nova surveys and more recent measurements of cluster abundances. 

Optical surveys show a correlation between SN la properties and the morphology of their 
host galaxy. While a population of brighter SN la with a slow luminosity decline is more 
common in late-type spiral and irregular galaxies with recent star formation (indicating a 
short delay time between their formation and the explosion), a faint er and more r apidly 
decaying populatio n of SN la is more common in early-typ e galaxies (Hamuv et al. 1996; 
Ivan ov et al.ll2000l) . This is clearly seen in the Fig. [35] from Gallaghe r et al. (2005) which 
shows the morphological type of the parent galaxy versus the decline rate of SNe la. SNe la 
in galaxies whose populations have a characteristic age greater than 5 Gyr are ~1 mag fainter 
than those found in galaxies with younger populations. Based on the evolution of the SN la 
rate with redshift, its c orrelation with the colour and radio luminosity of the parent galaxies, 
Mannucci et al. Ulooi) concluded that the delay time distribution of SN la is bimodal. About 
half of the SN la explode ~ 10 s years after the formation of the progenitor binary system and 
the delay time distribution of the other half of SN la could be described by an exponential 
function with a decay time of ~3 Gyr. This diversity should be reflected in the abundance 
yields. The brighter supernovae are expected to produce more Ni and less a elements (Si, S, 
Ar, Ca) than the fainter ones. 

Fig. [36] shows the logio( Fe/Si) against the logio( Fe/Q) for differe nt spatial regions i n 
clusters of galax ies Hydra A dSimionescu et al.ll2009h. Sersic 159-0 3 |de Plaa et al. l2006h. 
Centaurus (Matsushita et al. 2007a), Fornax (Matsushita et al. 2007b), M 87 {Matsushit a et al.l 
12003b . and A1060 ([Sato et al.ll2007r) . The plot shows a relatively large spread in the abun- 
dance patterns determined for different clusters. In particular the differences in the Fe/Si 
ratio between Hydra A and M 87 are striking. Each overplotted line in the figure sho ws a 
mixture of yields of a different SN la model with SNcc yields bv lTsuiimoto et all d 19951) in- 
tegrated over the Salpeter IMF. Using SNcc yields by different authors, assuming different 
progenitor pre- enrichment levels, or integrating the SNcc yields over a different IMF would 
move these curves (see Fig. 11 in lSimionescu et al1l2009l) . But an intrinsic variation in the 
metallicity or IMF of the SNcc progenitors alone could not explain the observed variety 
of abundance patterns. Thus it seems that, indeed, a variety of SN la flavors with different 
densities at the transition from the subsonic to supersonic flame velocities is needed. 

While the ratio of the a elements to Fe can be different for different SNe la, the ratio 
between the a elements is very similar for all explosion mechanisms, because the fuel (C/O) 
is always the same and the burning conditions where the a elements are synthesized do 
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Fig. 36 Intrinsic variation of cluster abundances and their average yields. Left panel: The observed [Fe/Si] 
versus [Fe/O] for different spatial regions of a sample of six clusters of galax ies. The curves indicat e mod- 
els of a mixture of SNcc prod ucts assuming a Salpe ter initial mass function iTsujimoto et al. 1995) mixed 
with different SN la products ilwamoto et al. 1999). The dotted line connects points of the same number 
contributions of SN la to the enrichment of th e ICM. The abundance patterns in the Hydra cluster and in 
M 87 favor very different SN la models. (After Simionescu et al. 2009). Right panel: Reconstruction of the 
avera ge abundance patter ns observed in a sample of 22 clusters of galaxies with theoretical supernova yields 
(from de Plaa et al. 2007j). T he model describing th e data the best is a combination of SN la yields from a ID 
de layed detonation mod el bv lBadenes et alj )2006h and SNcc yields for progenitor s with Solar metallicities 
bv lNomoto et alj 120061) integrated over the Salpeter IMF (from lde Plaa et all2007l) . 



not vary much in SN la. Therefore the relative abundances of the a elements, especia lly the 
Ca/Ar ratio, are the best diagnostic of the quality of SN la models, de Pl aa et all d2007l) anal- 
ysed the abundance patterns of Si, S, Ar, Ca, Fe, and Ni in a sample of 22 clusters observed 
with XMM-Newton within Q.IRsqq. Because the SNcc do not have a significant impact on 
the Ar/Ca ratio, they used the abundances of these elements to find the best fitting SN la 
model. Delayed detonation models clearly provide a better fit to the data tha n the deflagra- 
tion scenarios, but none of the supe rnova models by Ilwamoto et alj d 19991) fit the Ar and 
Ca abundances, de Pl aa et al. I d2007t) obtained a good fit wit h a one-d imensional delayed- 
detonation SN la model, calculated on a grid introduced in iBadenes et all d2003h . with a 
deflagration-to-detonation density of 2.2 x 10 7 g cm -3 and kinetic energy of 1. 16x 10 51 ergs, 
which was shown to fit best the properties of the Tycho supernova remnant jBadenes et alj 
l2006h . 

Here it is interesting to note that Tycho seems to be a remnant of a fairly typical SN la 
with an average luminosity and the data points in the left panel of Fig.[36]also tend to be cen- 
tred a round the WDD2 model which is the closest to the best fitting model of Ba denes et alj 
d2006h . Several other observa tional clues point to ward the delayed deto nation models. Opti- 
cal light curves (Woosley et al. 2007) a nd spe ctra (Mazzali et al. 2007), the X-ray emission 
of SN la remnants (Bade nes et alj2006l ,l2008) all show that the structure of the SN la ejecta 
has to be very close to what one-dimensional delayed detonation models predict: most of 
the white dwarf is burnt, Fe-peak nuclei stay in the inner ejecta and the intermediate mass 
elements stay in the outer ejecta. The variation of the peak brightness, which correlates with 
the production of 56 Ni and anti-correlates with the production of Si-group elements, can 
also be explained in the framework of the delayed detonation models by a variation of the 
deflagration-to-detonation transition density (transition from subsonic to supersonic flame 
velocities). 

As shown by Ide Plaa et al.1 fc007l) the relative abundances of chemical elements in the 
ICM can, in principle, be also used to distinguish between SNcc models with different level 
of pre-enrichment of the progenitors and with different initial-mass functions (IMF). Com- 
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paring the measured abunda nces with the SNcc yields predicted by the model of lNomoto et al.l 
(2006D. lde Plaa et aD J2007h found that the dominant fraction of SNcc progenitors contribut- 
ing to the enrichment of the ICM was already enriched and their mass distribution was close 
to the Salpeter IMF. The main constraints on the SNcc progenitors are based on the mea- 
sured O/Ne ratio in the ICM. 

While elements from O up to the Fe-group are primarily produced in supernovae, car- 
bon and nitrogen are believed to originate from a wide range of sources: winds of short-lived 
massive metal rich stars, l onger lived low- and intermediate-mass stars, and an early gener- 
ation of massive stars (e.g. iGustafsson et"allll999l ; IChiappini et ai]|2003l ; iMevnet & Maederl 
l2002h . The ICM abun dance of C and N was measured only in a few nearby bright elliptical 
galaxies and clusters ( IPeterson et al.ll200 3: Werner et al. l2006al) . The small O/Fe ratio and 
large C/Fe and N/Fe found in M 87 bv lWerner et all ( l2006ah suggest that the main sources 
of C and N are not the massive stars that also produce large quantities of O, but the low- and 
intermediate-mass asymptotic giant branch stars. 



5.4 Metallicity evolution and differences between clusters 

As we mentioned before, cluster data obtained by ASCA were consistent with no evolu- 
tion in the ICM metallicity out to redshift z ~ 0.4. More recent observations with Chan - 
dra, however, indica te a significant evolution tealestra et alj|2007l ; iMaughan et al.l l2008). 
Balestra et al. (2007), who investigated the Fe abundance evolution on a sample of 56 clus- 
ters, found that while in the redshift range of z — 0.3-0.5 the average IC M Fe abundance 
is ^ .4 Solar, above redshift z ~ 0.5 the metallicity drops to - 0.25 Solar. Maughan et al.l 
(2008) looked at a larger sample of 1 16 clusters at 0.1 < z < 1.3 in the Chandra archive and 
essentially confirmed the results of iBalestra et al.l ( 120071) . They found that the abundances 
drop by ~50% between z ~ 0. 1 and z ~ 1 and the evolution is still present if the cluster cores 
(the inner O.I5R500) are excluded from the analysis. This result indicates that the abundance 
drop is not due to the lack of strong cool cores at large redshifts. 

Balestra et"aD ( 120071) also found a trend of the Fe abundance with the cluster tempera- 
ture. Within (0.15-0.3)/? v ir in clusters below 5 keV the Fe abundance is on average a factor 
of ~2 larger than in the hotter clusters. The Fe abunda nce values m e asured within O.2.R500 
for the sample of low redshift 22 clusters analyzed by Ide Plaa et al. show the same 

trend with the cluster temperature: while in hot massive clusters (kT > 5 keV) the Fe abun- 
dance seems to be constant and equal to ~ 0.3 solar, for cooler clusters, in the temperature 
range of 2-4 keV, the Fe abundance shows a range of values between 0.2-0.9 solar. This 
trend is probably linked to the changing stellar mass over gas mass ratio in clusters. While 
the gas mass fraction was found to, in general, increase with the clust er temperature with 
a large spread of values especially below 5 keV dSanderson et al. | |2003|). the gas fr action in 



clusters above > 5 keV is constant with the ICM temperature I Allen et al„ 2008). The Fe 



mass ove r light ratio seems to be constant as a function of temperature (e.g. iRenzinill 19971) . 
iLin et al.l ( 120031) find that the stellar mass fraction decreases by a factor of 1.8 from low- to 
high- mass clusters, and the ICM to stellar mass ratio increases from 5.9 to 10.4. This trend 
suggests a decrease of star formation efficiency with increasing cluster mass and provides a 
natural explanation for the on average lower iron abundance in the more massive clusters. 
The large range of Fe abundance values in the clusters with lower temperatures might be 
due to the large spread in the gas mass fraction of the lower mass clusters and due to the 
range of enrichment times and gas densities in the cool cores of clusters. 
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Fukazawa et al.l Jl998h and lBaumgartner et all d2005h showed based on ASCA data an 
increase of the Si/Fe ratio from 0.7 to 3 Solar in the 2-8 keV temperature range. This result 
was not confirmed by later XMM-Newton observations. The abundance ratios of Si/Fe, S/Fe, 
Ca/Fe, and Ni/Fe measured within .2^00 in the sample of 22 clusters with temperatures in 
the range of 1.5-10 keV analyzed bv lde Plaa et"aD d2007h are consistent with being constant 
as a function of temperature. The intrinsic spread in the abundance ratios is smaller than 
30%. This result shows that the abundance patterns do not change with the observed drop of 
star formation efficiency with the cluster mass and the ratios of SN la and SNcc contributing 
to the ICM enrichment are similar in clusters of all masses^. 

The observed drop of the Fe abundance with increasing redshift could be at least partly 
connected to the metallicity-temperature relation. The clusters observed at high redshift are 
primarily massive clusters which eventually evolve into hot clusters like those with lower 
abundances at low redshifts. The evolution of the metallicity below redshift z ~ 1 therefore 
might be substantially weaker than the observations suggest. 

5.5 How good are the assumptions behind the abundance measurements? 

The chemical abundances in the ICM are determined assuming a certain plasma model, a 
certain ICM temperature structure, optically thin plasma, and solar abundances for elements 
lighter than O or N. How good are these assumptions and how can they bias our abundance 
determinations? 

5.5.7 How much bias due to spectral modeling? 

Usua l ly two plasma models are used for fitting cluster spectra: ME KAL ( lMeweetal.ll 19851 , 
Il986l : lKaastr3l 19921 : lliedahl et al.lll995h and APEC dSmith et alJboOlh . A more updated 
version of MEKAL com pared to the one a vailable in XSPEC is implemented in the SPEX 
spectral fitting package dKaastra et alJl99rjh . MEKAL and APEC treat differently the atomic 
physics and differences are also in their line libraries. But despite small differences in the 
best fit temperature structure, the abund ances derived with the two plasma codes are gen- 
erally consistent within the errors llBuoteetal .1120031 : 1 Sanders & FabiarJl2006l : lde Plaa et all 
l2007h . The two plasma codes thus do not show any systematic inconsistencies. The imple- 
mented ionization balance for Ni is inaccurate and can significantly bias the Ni abundance 
measurements; an update is planned shortly for SPEX (J.S. Kaastra, private communica- 
tion). 

One has to be careful and keep in mind possible biases when interpreting the results of 
measurements of chemical abundances in clusters. The most serious biases arise when the 
temperature structure in the extraction region is not modeled properly. Historically, X-ray 
observations showed sign ifica ntly subsolar a bundances in galaxies and groups of galaxies. 
iBuote & Canizaresl dl994l) and lBuotel d2000bl) showed that a low resolution spectrum which 
is composed of intrinsically two temperature components with similar emission measures, 
one below 1 keV and one above 1 keV, with an average temperature of about 1 keV fitted 
with a single temperature model will seriously bias the Fe abundance determination. Since 
the lower temperature component preferentially excites emission lines in the Fe-L complex 

2 The increase of the Si/Fe ratio with the cluster temperature observed by ASCA might have been due to a 
bias introduced by the strong temperature gradients in the hot clusters with cooling cores. If the temperature 
of the best fit plasma model (the average cluster temperature) is higher than the temperature of the cooler gas 
with high Si line emissivity, then the model will naturally predict a higher Si abundance than the actual value. 
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Fig. 37 Left panel: The influence of the temperature model on the best fit Fe abundance as seen in the radial 
Fe abundance distribution in NGC 1399 modeled with single- and two-temperature models (from Buote 
2002). Right panel: The helium abundance is usually assumed to be Solar in clusters. If the helium abundance 
in cluster cores is increased by sedimentation and we do not take it into account in the spectral modeling, our 
best fit abundances will be lower than the real values. The effect is indicated in the plot for simulated spectra 
of a 4 keV cluster. The simulation was made assuming a 100 ks observation with XMM-Newton EPIC and a 
circular extraction region with a radius of 1'. 



below 1 keV (Fe xii-xxi) and the higher temperature component excites the Fe-L lines in 
the range of ~ 1-1.4 keV (Fexxi-xxiv) the shape of the Fe-L complex will be flatter than 
that produced by a single-temperature model with the average temperature of the two com- 
ponents (see also Fig. 123) . The single-temperature model can fit the flat shape of the Fe-L 
complex only with a relatively low Fe abundance and higher bremsstrahlung continuum. 
The obtained Fe abundance is thus much lower than the actual value (see the left panel of 
Fig.O. 

In the past few years several authors (e.g. lBuote et alj2003MWerner et alj2006bl : lMatsushita et al.l 

2007a) showed that in the complex cool cores of groups and clusters the obtained abun- 
dances of essentially all elements increase when two-temperature models or more com- 
plicated differential emission measure (DEM) models are used instead of simple single- 
temperature models. 

For a multi-temperature plasma with average temperatures of ~2^1 keV, where both 
Fe-L and Fe-K emission lines are seen with s imilar statistics, fi t ted w ith a single temper- 
ature model an "inverse Fe bias" is observed. Simi onescu et al. point out that in a 
multi-temperature plasma at these temperatures the cool component produces stronger Fe-L 
lines and the hotter component emits more Fe-K emission than that expected for the single- 
temperature plasma with the average temperature. Therefore, by fitting a single temperature 
model to such p lasma one will obtain a higher Fe abundance than the actual value (see also 
iRasia et alj2008h . 

These results show that for chemical abundance studies in the ICM a good description 
of the temperature structure is crucial. Accurate abundance determinations require deep ob- 
servations which allow to test different differential emission measure distributions to obtain 
the best description of the true temperature structure in the cluster cores. 
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5.5.2 Possible biases due to helium sedimentation? 

In the gra vitational potential of clusters helium and other metals can diffuse toward the clus - 
ter center (Fabian & Pring lej l 19771 : lAbramopoul os et al ]ll98ll : iGilfanov & Svunvaevll 19841) . 
The sedimentation is more efficient for the lighter elements and in a H-He plasma, the he- 
lium is expected to be more centrally peaked than hydrogen. The helium abundance can not 
be directly measured. The cluster spectra are commonly modeled under the assumption of 
solar or primordial He abundance values in the ICM. In the inner regions of galaxy clus- 
ters where the effects of sedimentation are the most important, the underestimation of the 
am ount of helium c an cau se inaccurate determination of gas densities and metallicities. 

IPeng & N agai (2009) solved a set of flow equatio ns for a H-He plasma for clusters 
assuming temperature profiles by Vikhl inin et al.l |2006). They showed that at 0.06r 5 oo of a 
1 1 Gyr massive (10 keV) cluster the He/H abundance ratio can be as high as 4.3 Solar. These 
predictions are, however, made for unmagnetized plasma assuming no turbulent suppression 
of par ticle diffusion and no m ixing in the ICM, therefore they should be treated as upper 
limits. lEttori & Fabianl foood) pointed out that for an actual helium abundance of ~4 times 
the value assumed in the spectral modeling, would cause an under-determination of the Fe 
abundance by a factor of 2. As we show in the righ panel of Fig.[37]the same is approximately 
true for the other elements as well. 

5.5.3 Resonance scattering in cooling cores? 

For th e strongest resonant lines the ICM can be moderately optically thick dGilfanov et al] 
Il987h and as we already mentioned in Sect. 3.5 resonant scattering at the energies of these 
lines changes the observed line intensities. This can affect the measurements of elemental 
abu ndances in the dense core s of clusters. 

ISanders & Fabianl ( Eooih constructed a model that takes into account resonant scattering 
for thousands of resonance lines. They tested their model on Chandra data of two clusters 
with central abundance dips. They conclude that metallicities in cluster cores could be un- 
derestimated due to resonant scattering by at most 10%. Resonance scattering can therefore 
not be the cause for the observed central abundance dips. 

In order to determine the amount of resonant scattering I Werner et aD|2009) studied the 
line ratios of the Fexvn ion in 5 elliptical galaxies (see Sect. 3.5). They concluded that 
neglecting the effects of resonant scattering in the spectral fitting of the inner 2 kpc core of 
NGC 4636 will lead to underestimates of the chemical abundances of Fe and O by ~ 10- 
20%. 



6 Beyond Thermal Equilibrium 

The intra-cluster medium is usually assumed to be in collisional ionization equilibrium 
(CIE). Plasma in CIE is optically thin for its own radiation and external radiation fields 
do not affect the ionization balance which is entirely determined by the temperature of the 
plasma. The ionization and recombination rates in CIE plasma have come to a balance and 
the electron and ion temperatures are in equilibrium. The energy distribution of the electrons 
in the ICM is usually described by a Maxwellian. 

These assumptions are mostly justified but there are certain situations when the plasma 
deviates from the thermal equilibrium. The shock heated plasma in the cluster outskirts with 
densities n e < 10~ 4 cm -3 is likely to be out of ionization balance. The equilibration time 
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scale at those low densities is longer than the age of the cluster and therefore the gas may 
still ionize. At the lowest density parts of warm-hot intergalactic medium filaments pho- 
toionization by the diffuse radiation field of galaxies and the cosmic background becomes 
important. Deviations from a Maxwellian electron distribution could occur in the cooling 
cores of clusters where multiple gas phases are present. If electrons diffuse from a hotter 
into a cooler phase they over-ionize the gas. Plasma could also be out of equilibrium in the 
vicinity of shocks, which may alter the electron/ion equilibrium. Shocks could also accel- 
erate electrons and produce non-Maxwellian tails in the ICM electron distribution. Future 
observatories with large photon collecting area and high spectral resolution will be able to 
study the spectral signature of these non-equilibrium effects in detail. Moreover, radio ob- 
servations of clusters reveal the presence of relativistic electrons emitting synchrotron radio 
emission. This relativistic plasma in clusters may also emit Inverse-Compton radiation in 
the form of a power law shaped X-ray continuum. 

Below we briefly describe what are the spectral signatures of non-equilibrium effects in 
the multiphase plasma and around shocks, and we summarize the state of observations of 
X-ray emission from relativistic plasma. 



6.1 Multiphase plasma 

High signal-to-noise X-ray spectra re veal that multiple spatially unresolved gas phases are 

presen t in the cooling cores of clusters dKaastra et alj2004llSanders & Fabiarj2007l ; lsimionescu et alj 
l2008bl) . In M87 and Perseus, the most nearby and the bri ghtest cooling cores , images ob- 
tained by Chandra reveal filaments and blobs in soft X-rays ( Forman et al .l2007l : lFabian et alJ 
2006). But the multi-phase extends to much lower temperatures. Narrow band optical ob- 
servations of cooling cores reveal Ha emitting warm (10 4 K) gas which both in M87 and i n 
Perseus is often associated with soft X-ray filaments dSparks et alj2004l ; lFabian et alj2006h . 
The Ha filaments i n Perseus are composed o f gas with temperatures of 10 4 K down to 50 K 
dHatch et"al]|2005h . lSanders & Fabiarj d2007h showed that they are surrounded and possibly 
mixed with soft X-ray emitting plasma with temperatures of 0.5-1 keV. 

If electrons penetrate from the hotter to the cooler phase they over-ionize the "cool" 
plasma which may result in observable low energy X-ray line emission. In Fig.[38]we show 
a simulated EPIC/pn spectrum (100 ks) of the Eastern X-ray arm of M87. Within this ex- 
traction region, we assume the presence of a filament with a diameter of 60 pc and length 
of 3.5 kpc filled with 10 4 K gas with an average density of n e = 10 cm -3 . We assume 
that 0.05% of all the electrons within the filament have penetrated there from the ambient 
ICM and have a Maxwellian energy distribution with a temperature of 2 keV. As the spec- 
trum in Fig. [38] shows, a deep observation with a well calibrated instrument would even at 
medium spectral resolution reveal the line emission from Ne m-vi at 0.85-0.9 keV, O iv-vn at 
-0.55 keV, Nm-vi between 0.40-0.42 keV, Civ-v at 0.35, 0.37 keV, and at 0.30-0.31 keV. 
This exercise demonstrates the potential of X-ray spectroscopy in constraining the diffusion 
of electrons from the hot ICM phase into the warm filaments. In reality such attempts are 
hampered by the complicated multi-temperature structure of the ICM around the filaments. 
Detailed studies of the electron diffusion will be possible with future micro-calorimeters or 
transition-edge sensors sensitive at low energies. 
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Fig. 38 Simulated EPIC/pn spectrum of the Eastern X-ray arm of M87 assuming a presence of a filament 
with a diameter of 60 pc and length of 3.5 kpc filled with 10 4 K gas with an average density of n e = 10 cm -3 . 
We assumed that 0.05% of all the electrons within the filament penetrated there from the ambient ICM and 
have a Maxwellian distribution with a temperature of 2 keV. The model plotted with the red line shows 
the emission from the ambient hot ICM. Assuming a well calibrated instrument the excess line emission at 
~0.85 keV would be clearly revealed. 



6.2 Shocks and non-equilibrium 

In the downstream regions of shock fronts the electron and ion temperatures of the plasma 
might be different and the ICM is out of ionization equilibrium. 

In a collisional plasma, protons are heated dissipatively at the shock layer which has 
approximately the width of the collisional mean free path, but the faster moving electrons do 
not feel the shock (for Mach numbers smaller than ~ 40), they get adiabatically compressed 
and equilibrate with the ions via Coulomb collisions reaching the post-shock temperature 
predicted by the Rankine-Hugoniot jump condition. Astrophysical shocks in a magnetized 
plasma, like the ICM, are expected to be "collisionless". Which means that the electron and 
ion temperature jump occurs on a spatial scale much smaller than the collisional mean free 
path, but ions and electrons are heated by the shock to different temperatures. In heliospheric 
shocks with moderate Mach numbers electro ns are heated much l ess than protons, barely 
above the adiabatic compression temperature dSchwartz et al.|[l988h . The plasma behind the 
shocks is thus expected to be in electron/ion temperature non-equilibrium. Assuming that 
no kinetic energy goes into acceleration of cosmic rays, the post-shock temperature after the 
Coulomb equilibration time scale reaches the same temperature as the one predicted by the 
Rankine-Hugoniot jump condition. 

iMarkevitchl d2006h used the bow shock in the "Bullet cluster" (IE 0657-56) to deter- 
mine whether the electrons in the intra-cluster plasma are heated directly at the shock to the 
equilibrium temperature (instant-equilibration), or whether they equilibrate via Coulomb 
collisions with protons. They measure the temperature in front of the shock front seen as a 
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Fig. 39 Left panel: Electron-ion equilibration in IE 0657-56 (the "Bullet cluster"). The two data points show 
the deprojected electron temperature for two narrow post shock regions. The data points are overlaid on the 
model predictions for instant equilibration (light gray) and for adiabatic compression followed by collisional 
equili bration (dark gray) . The velocity of the post-shock gas relative to the shock front is ~ 1600 km s . 
(From Markevitch 2006). Right panel: The 6.9-7.0 keV part (rest frame energies) of a simulated Astro-H 
micro-calorimeter spectrum of the downstream region of a M = 2.2 shock propagating through a 10 8 K 
plasma. The Fe xxvi satellite line is clearly stronger (with a significance of 3<r for these simulations) than 
that predicted by th e thermal model with a Maxwellian electron distribution indicated by the full line. (From 
lKaastraetaIT2 009). 



gas density jump, from which, using the adiabatic and Rankine-Hugoniot jump conditions, 
they determine the expected post shock adiabatic and instant-equilibration electron temper- 
atures which they directly compare with observed data. The shock in the "Bullet cluster" is 
happening approximately in the plane of the sky, which means that the downstream velocity 
of the shocked gas spreads out the time dependence of the electron temperature in the plane 
of the sky. The measured temperatures behind the shock in the bullet cluster are consistent 
with instant heating of electrons. Adiabatic compression followed by equilibration on the 
collisional time scale is excluded at the 95% confidence level (see Fig. 139b- 

As the temperature of the plasma suddenly rises at a shock front, it gets out of ionization 
balance. While the ionization state of ions still reflects the pre-shock temperature of the 
plasma after the instant heating, the electron temperature is higher. The plasma is under- 
ionized compared to the equilibrium case and the ionization balance must be recovered 
by collisions. Until the ionization balance recovers there will be more ionization in the 
plasma than recombination. The exact ionization state of shocked plasma in non-equilibrium 
conditions depends on its temperature, density, and the time since it has been shocked. Such 
situation is often seen in supernova remnants. 

Shocks and turbulence in the magnetized ICM are believed to be sid es of non-thermal 
particle acceleration in cl usters (for reviews on particle acceleration see e.g. lPetrosian & B vkov 
2008; Bvkov et al. 2008). Stochastic t urbulence in the ICM ca n create a non-thermal tail to 
the Maxwellian electron distribution (Bvk ov &Uvarovll 19991) . The low energy end of the 
power-law electron distri bution enhances the i onization rates and modifies the degree of the 
ionization of the plasma. IPorquet et al.l (200 1 ) showed that the plasma is always more ion- 
ized for hybrid (Maxwellian plus power-law) electron distribution than for a Maxwellian 
distribution and the mean charge of a given element at a given temperature is increased. The 
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effect is more pronounced at lower temperatures. In groups of galaxies, low mass clusters or 
in cooling cores a peculiar ionization state of Fe which is inconsistent with the temperature 
determined from continuum emission c an potentially be a g ood tool to reveal such non- 
thermal tails in the electron distribution. IKaastra et alj (2009) show that a non-Maxwellian 
tail in the electron distribution behind a shock front propagating through 10 s K plasma 
cannot be revealed by current detectors. They simulate the electron distribution and the re- 
sulting X-ray spectra of the downstream region of a M = 2.2 shock in a X-ray bright hot 
cluster and show that a good indicator of hard non-thermal electrons is the enhancement 
of the equivalent widths of satellite lines which may be possible to detect with the X-ray 
micro-calorimeters on the future Astro-H satellite or on the proposed International X-ray 
Observatory (IXO). In the right panel of Fig.[39j we show the 6.9-7.0 keV part (rest frame 
energies) of a simulated 100 ks Astro-H spectrum of the downstream region of a M = 2.2 
shock propagating through the ICM in a cluster similar to Abell 2029. The Fe xxvi Lya 
lines and the Fe xxvj satellite line are visible in the spectrum and the satellite line is clearly 
stronger than that predicted by the thermal model with a Maxwellian electron distribution 
indicated by the full line. 



6.3 Non-thermal X-ray emission from relativistic plasma 

Direct evidence for the presence of electrons accelerated up to relativistic energies comes 
from observations of radio emission in clusters. Large-scale diffuse extended radio emis- 
sion in the form of hal os or relics has been observed in about 50 known clusters of galaxies 
jFeretti & Giov annini 2007J). This radiation is associated with the ICM and has no connec- 
tion to the cluster galaxies. It is clearly produced by synchrotron emission by a population 
of relativistic electrons in the ICM. Since the energy density of the cosmic microwave back- 
ground (CMB) radiation (~ 4x 10" 13 erg cm -3 ) is higher than the energy density in the intra- 
cluster magnetic field (3 x 10" 14 (5/^G) 2 erg cm -3 ), these relativistic electrons will radiate 
away most of their energy via inverse-Compton scattering of th e CMB photons, producing 
a power-law shaped X-ray continuum emission (Rep haelilll977h . The detection of both the 
non-thermal X-ray emission and the radio emission would be a powerful tool, which would 
allow us to determine both the volume averaged magnetic field in clusters and the energy in 
the population of relativistic electrons. But for the observed radio fluxes, detectable hard X- 
ray emission can only be produced if the magnetic fields are of the order of 0. 1 /iG. Faraday 
rotation measurements indicate that the intra-cluster magnetic fields are much stronger - of 
the order of 1-10 /jG. The volume averaged magnetic field might, however, be weaker than 
the strong magnetic field measured by Faraday rotation along our line of sight. 

Non-thermal X-ray emission could in principle also originate in non-thermal brems- 
strahlung of shock accelerated supra-thermal electrons i n the ICM. How ever, such non- 
thermal bremsstrahlung phase must be very short lived. IPetrosiar] ( l200ll) points out that 
in order to explain the reported hard X-ray luminosity in the Coma cluster (see later) by 
bremsstrahlung emission, the continuous input of energy into the ICM would increase the 
ICM temperature to 10 10 K in Hubble time. Also the high-energy electrons > 50 keV can not 
be confined by the gravitational potential of the cluster and will escape in a crossing time of 
< 1.5 x 10 7 years, unless they are confined by the intra-cluster magnetic field. Alternatively, 
hard X-ray emission with power-law spectra of r ~ 1.5 could be produced by synchrotron 
emission from ultra-relativisti c electrons and p ositrons produced by the interaction of rela- 
tivistic protons with the CMB dlnoue et al . 2005). 
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The nearest hot massiv e cluster of galaxies with a big radio halo is the Coma cluster 
dFeretti & Giovaniunill 19981) . Therefore, Coma has been the primary t arget to identify non- 
therm al X-rays. T heir detection has been rep orted with both BeppoSAX dFusco-Femiano et al.l 
Il999h and RXTE (Rephaeli & Gruber 2003 ). Mor e recently, Coma has also been observed 
with the INTEGRAL satellite fackert et al.ll2007l) . A combined analysis of XMM-Newton 
and INTEGRAL data revealed the presence of hotter (~12 keV) gas in the south -west region 
overl apping with a radio halo. This hot gas was probab ly heated in a merger feckert et al.l 
|2007|) . Analysing INTEGRAL, ROSAT, and RXTE data, iLutovinov etall d2008h found that 
the global Coma spectrum is w ell fitted with a the rmal model and the evidence for a hard 
excess is very marginal (T.6tr). lAjello et aD J2009h analysed combined XMM-Newton and 
hard X-ray Swift-BAT spectra and even though they could not rule out the presence of non- 
thermal hard X-ray emission outside of their 10' region, within their field of view they found 
a good fit with two thermal models, with no need for a non-thermal component. Iwik et al.l 
(2009) analyzed combined Suzaku Hard X-ray Detector (HXD-PIN) and XMM-Newton mo- 
saic spectra of a larger 34'x34' region centered on the Coma cluster. They found no statis- 
tically significant evidence for non-thermal emission implying a lower limit of 0.15 yuG on 
the cluster averaged magnetic field. 

Suzaku observations of the bright radio relic in the merging cluster Abell 3667 show 
that, at least in some cases, magnetic fields can be high even at distances of ~ 1 Mpc from 
the cluster core. The upper limits on the non-thermal emission from Suzaku put a lower 
limit of 2.3 yuG on the magnetic field in the radio relic. The non-th ermal energy density i n 
the relics is >7% of th e thermal energy density and likely near 20% (Nakazawa et al. 2009). 

lAiello etaT] ( l2009h analysed the spectra of 10 clusters of galaxies detected above 15 keV 
with the Swift-BAT. Except of the Perseus cluster the spectrum of which is probably con- 
taminated by hard X-ray emission from the central AGN, the spectra of all clusters are well 
fitted by a simple thermal model. Their stacked spectrum of 8 clusters (except Perseus and 
Coma) also confirms the absence of any non-thermal high energy component down to a flux 
of 1.9 x 10" 12 erg cirrV 1 in the 50-100 keV band. 

Although there have been several other reports claiming the detection of both hard and 
soft X-ray emission in excess to the ICM emission , interpreted as bein g possibly of non- 
thermal origin (for reviews see Rephaeli et al. 200jl; iDurret et al.ll2008l) . a solid identifica- 
tion of non-thermal X-rays from clusters is still lacking. In the relatively near future, X-ray 
satellites with imaging hard X-ray optics sensitive up to 40 keV will be launched: Astro-H, 
NuSTAR, Simbol-X. They will shed more light on the presence of very hot thermal gas, non- 
thermal and supra-thermal electron populations, and on intra-cluster magnetic fields. The 
Fermi Gamma- ray Large Array Space Telescope (GLAST) might detect gamma ray emis- 
sion associated with relativistic intra-cluster ions. These missions promise a big progress in 
our understanding of the non-thermal particle population in clusters of galaxies. 



7 Future Potential of X-ray Spectroscopy 

One of the most awaited features of the future of X-ray observatories is high resolution 
spectroscopy of the next generation of imaging devices based on calorimetric measurements 
of the photon energy in detectors cooled to few tens of mK. The two missions for which such 
detectors are envisioned and actually play a major role in the scope of these missions are 
Astro-H in Japan (with a possible launch in 2013^1) and the International X-ray Observatory 



3 http://astro-h.isas.jaxa.jp/ 
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Fig. 40 Simulated 1X0 NFI spectrum of cluster emission with a flux of 3 ■ 10~ 13 erg s - ' cm -2 in the NFI 
field-of-view at a redshift z ~ 0.2 with an exposure of 100 ksec. The inset shows the spectral structure of 
major Fe L-shell lines around 1 keV. Spectral line broadening for plasma with mean turbulent velocities of 
100 km/s and 500 km/s are compared in this inset. 

{IXO) currently discussed by ESA, NA SA, and JAXA as a mission for about 2018-20210 
While Astro-H dTakahashi et al.ll2009l) is already being built, the IXO mission has not yet 
been fully approved and is still in the first phase of technical studies. Nevertheless, to give an 
outlook on the coming spectroscopic capabilities of these future X-ray missions, we will use 
the current prospective response files for the IXO mission with the planned so called narrow 
field detectors (NFI), which provide high resolution spectroscopy of about 2-3 eV energy 
resolution, to illustrate the prospects. In addition to providing high resolution spectroscopy 
IXO will also significantly increase the photon collecting power compared to XMM-Newton 
with a proposed collecting area around 1 keV of about 3 m 2 to provide a better photon 
statistics to fully exploit the higher spectral resolution statistically. It will also increase the 
redshift range of cluster studies providing a detailed picture of galaxy cluster astrophysics 
up to redshifts of 2. 

With the prospective energy resolution it will be possible to trace gas motions in the ICM 
with an accuracy of better than 100 km s in the Fe lines in the X-ray spectra. This is made 
possible due to the fact that the thermal line broadening for lines from Fe ions is smaller by 
a factor of about 7 than the sound velocity which corresponds closely to the thermal velocity 
of the bulk of the baryons in the ICM, the protons. Fig.|40]shows the simulated spectrum as 
observed with the IXO NFI of a patch in a cluster with an aperture radius of 1 arcmin and a 
flux of 3 x 10" 13 erg s" 1 cm" 2 (at 0.5-2 keV). With an exposure of about 100 ks (40 ks) a 
Gaussian shape turbulent velocity broadening can be measured with an accuracy of 40 (70) 
km s _1 . Even for a cluster at redshift 1 with a flux of 10" 14 erg s -1 cm -2 velocity broadening 
accuracy would still be 100 km s _1 in an exposure of 100 ks. This will provide many new 
insights into the ICM physics, answering questions such as for example: how much turbulent 
energy is injected into the ICM by the interaction of cool cores with the central AGN? How 
much of the cool core heating comes from the dissipation of this turbulent energy? How is 



4 http://ixo.gsfc.nasa.gov/ 
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Fig. 41 Simulated IXO NFI spectrum of a galaxy group with an ICM temperature of 2 keV at redshift 2. Also 
shown is the estimated background spectrum. ICM properties like the temperature and abundance measure- 
ments can well be determined from this simulated exposure of 250 ksec. See text for details. 



the energy dissipated in merging clusters and how is this reflected in the turbulent and bulk 
motion structure of the ICM? How much can the mass estimate of clusters be improved if 
we can measure the energy content of the ICM in turbulent energy? 

Finally Fig. [41] illustrates the potential of IXO to provide detailed information on very 
distant clusters. Common objects at redshift 2 having a sky density of about 1 per deg 2 are 
poor clusters with masses of about 3 x 10 13 M Q and probably temperatures around 2 keV. 
Fig. 14 II shows t he a spectrum of su ch an object observed with the IXO NFI instrument for 
250 ks (see also lArnaud et al.ll2009t) . It will allow the measurement of the bulk temperature 
with an accuracy better than 5% and provide abundance determinations with interesting 
precision: Fe (±11%), Si (±18%), O, Mg (±30%). This will allow us to study the thermal and 
chemical enrichment history of galaxy clusters over most of the interesting cosmic epochs in 
which clusters existed. For the more nearby bright clusters IXO will very accurately measure 
the abundances of many chemical elements, including trace elements like Cr, Ti, Mn, or Co. 
The other elements, which we detect in clusters with XMM-Newton and Chandra (O, Ne, 
Mg, Si, S, Ar, Ca, Fe, Ni) will be traced with an unprecedented accuracy, which will help 
us to put better constrains on supernova yields and thus on their explosion mechanisms. 
Furthermore, the wide field imager on IXO will allow us to resolve the 2-dimensional metal 
distribution down to the relevant mixing physical scales and study in detail the process of 
metal injection. 

The filamentary warm-hot intergalactic medium (WHIM) permeating the cosmic web 
might contain up to 50% of baryons at redshifts of z < 2. This intergalactic gas is heated 
to temperatures between 10 5 and 10 7 K as it gathers in regions of overdensity 10-100 (with 
respect to the mean baryon density of the Universe; at the current epoch the mean baryon 
density is (ny) = 2 x 10" 7 cm -3 ). Because this gas has such a low density and because it 
is so highly ionized, it is very hard to detect. The current instrumentation only allows us 
to probe the highest overde nsity regions of the WHIM, such as the filament between the 
clusters Abell 222 and 223 dWerner et al1l200ot) . or to look for absorption features toward 
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bright continuum sources i n sight lines with k nown large scale structure, such as H2356-309 
behind the Sculptor Wall dBuote et al]|2009n . The combination of the large collecting area 
of IXO and the high spectral resolution of its diffraction grating spectrometers will allow us 
to begin to probe the low-overdensity regions of the intergalactic medium. High-resolution 
X-ray absorption spectroscopy of sight lines towards bright continuum sources may allow 
us to probe overdensities of 10 or more (assuming a metallicity of 0.1 Solar), through the 
absorption features of H- and He-like ions of oxygen. Emission line imaging spectroscopy 
of the WHIM should be pursued with a large field-of-view (of order one degree), short focal 
length, large grasp (effective area times field-of-view) soft X-ray telescope, with a high 
resolution imaging X-ray spectrometer at the focus. A successful experiment will need a 
grasp of order sever al hundred cm 2 deg 2 , and a spectral resolution of order 1-2 eV or better 
JPaerels et al1l2008h . 

Apart from these expectations, a new mission like IXO capable of opening a new param- 
eter space for X-ray spectroscopy and simultaneously increasing the sensitivity substantially 
and thus the reach of the observations, has an enormous potential for new discoveries. The 
most important discoveries of such missions generally come as surprises. Since the concept 
of IXO is already realistic - apart from detailed technical problems still to be solved - and 
also not much more costly than previous missions, we can expect its realization. Therefore 
we can foresee a bright future for X-ray spectroscopic research also over the coming two 
decades, and clusters will surely be among the most interesting objects to study. 

Acknowledgements We like to thank A. Bykov, G. Pratt, Carles Badenes, and Aurora Simionescu for a 
critical reading of the manuscript and for discussions. NW acknowledges support provided by the National 
Aeronautics and Space Administration through Chandra Postdoctoral Fellowship Award Number PF8-90056 
issued by the Chandra X-ray Observatory Center, which is operated by the Smithsonian Astrophysical Obser- 
vatory for and on behalf of the National Aeronautics and Space Administration under contract NAS8-03060. 
HB acknowledges support for the research group through The Cluster of Excellence "Origin and Structure 
of the Universe", funded by the Excellence Initiative of the Federal Government of Germany, EXC project 
number 153 and support from the DfG Transregio Programme TR33. 



References 

Abramopoulos, R, Chanan, G. A., & Ku, W. H.-M. 1981, ApJ, 248, 429 
Ajello, M., Rebusco, P., Cappelluti, N., et al. 2009, ApJ, 690, 367 
Akritas, M. G. & Bershady, M. A. 1996, ApJ, 470, 706 
Allen, S. W. & Fabian, A. C. 1998, MNRAS, 297, L63 

Allen, S. W., Fabian, A. C, Johnstone, R. M., Arnaud, K. A., & Nulsen, P. E. J. 2001, 
MNRAS, 322, 589 

Allen, S. W., Rapetti, D. A., Schmidt, R. W., et al. 2008, MNRAS, 383, 879 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Arnaud, M., Barcons, X., Barret, D., et al. 2009, Experimental Astronomy, 23, 139 

Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2005, A&A, 441, 893 

Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2007, A&A, 474, L37 

Arnaud, M. & Rothenflug, R. 1985, A&AS, 60, 425 

Arnaud, M., Rothenflug, R., Boulade, O., Vigroux, L., & Vangioni-Flam, E. 1992, A&A, 
254,49 

Badenes, C, Borkowski, K. J., Hughes, J. P., Hwang, U., & Bravo, E. 2006, ApJ, 645, 1373 
Badenes, C., Bravo, E., Borkowski, K. J., & Dommguez, I. 2003, ApJ, 593, 358 
Badenes, C., Hughes, J. P., Cassam-Chenai, G., & Bravo, E. 2008, ApJ, 680, 1 149 



57 



Balestra, I., Tozzi, P., Ettori, S., et al. 2007, A&A, 462, 429 

Baumgartner, W. H., Loewenstein, M., Horner, D. J., & Mushotzky, R. F. 2005, ApJ, 620, 
680 

Belsole, E., Sauvageot, J.-L., Pratt, G. W., & Bourdin, H. 2005, A&A, 430, 385 
Birkinshaw, M. 1999, Phys. Rep., 310, 97 

Bohringer, H. 1998, in The Local Bubble and Beyond, D.Breitschwerdt, M.J. Freiberg, J. 

Triimper (eds.), Springer Verlag, Heidelberg, 341 
Bohringer, H., Belsole, E., Kennea, J., et al. 2001, A&A, 365, L181 
Bohringer, H. & Hensler, G. 1989, A&A, 215, 147 

Bohringer, H., Matsushita, K., Churazov, E., Finoguenov, A., & Ikebe, Y. 2004, A&A, 416, 
L21 

Bohringer, H., Matsushita, K., Churazov, E., Ikebe, Y., & Chen, Y. 2002, A&A, 382, 804 
Bohringer, H., Nulsen, P. E. J., Braun, R., & Fabian, A. C. 1995, MNRAS, 274, L67 
Borgani, S., Finoguenov, A., Kay, S. T., et al. 2005, MNRAS, 361, 233 
Borgani, S., Murante, G, Springel, V., et al. 2004, MNRAS, 348, 1078 
Borgani, S., Rosati, P., Tozzi, P., et al. 2001, ApJ, 561, 13 

Bourdin, H., Sauvageot, J.-L., Slezak, E., Bijaoui, A., & Teyssier, R. 2004, A&A, 414, 429 

Braglia, F, Pierini, D., & Bohringer, H. 2007, A&A, 470, 425 

Brickhouse, N. S., Raymond, J. C, & Smith, B. W. 1995, ApJS, 97, 551 

Briel, U. G. & Henry, J. P. 1994, Nature, 372, 439 

Briel, U. G, Henry, J. P., & Boehringer, H. 1992, A&A, 259, L31 

Buote, D. A. 2000a, ApJ, 539, 172 

Buote, D. A. 2000b, MNRAS, 311, 176 

Buote, D. A. 2002, ApJ, 574, LI 35 

Buote, D. A. & Canizares, C. R. 1994, ApJ, 427, 86 

Buote, D. A. & Fabian, A. C. 1998, MNRAS, 296, 977 

Buote, D. A., Gastaldello, F, Humphrey, P. J., et al. 2007, ApJ, 664, 123 

Buote, D. A., Lewis, A. D., Brighenti, F, & Mathews, W. G. 2003, ApJ, 595, 151 

Buote, D. A., Zappacosta, L., Fang, T., et al. 2009, ArXiv e-prints 

Bykov, A. M., Dolag, K., & Durret, F. 2008, Space Science Reviews, 134, 119 

Bykov, A. M. & Uvarov, Y. A. 1999, JETP, 88, 465 

Canizares, C. R., Clark, G. W., Jernigan, J. G, & Markert, T. H. 1982, ApJ, 262, 33 
Canizares, C. R., Clark, G. W., Markert, T. H, et al. 1979, ApJ, 234, L33 
Cappellari, M. & Copin, Y. 2003, MNRAS, 342, 345 

Cavagnolo, K. W., Donahue, M., Voit, G. M., & Sun, M. 2009, ApJS, 182, 12 
Chiappini, C, Matteucci, F, & Meynet, G. 2003, A&A, 410, 257 

Churazov, E., Bruggen, M., Kaiser, C. R., Bohringer, H., & Forman, W. 2001, ApJ, 554, 
261 

Churazov, E., Forman, W., Jones, C, Sunyaev, R., & Bohringer, H. 2004, MNRAS, 347, 29 

Conselice, C. J., Gallagher, III, J. S., & Wyse, R. F. G. 2001, AJ, 122, 2281 

Cox, D. P. 1972, ApJ, 178, 143 

Cox, D. P. & Daltabuit, E. 1971, ApJ, 167, 113 

Cox, D. P. & Tucker, W. H. 1969, ApJ, 157, 1157 

David, L. P., Jones, C, & Forman, W. 1996, ApJ, 473, 692 

David, L. P., Nulsen, P. E. J., McNamara, B. R., et al. 2001, ApJ, 557, 546 

De Grandi, S., Ettori, S., Longhetti, M., & Molendi, S. 2004, A&A, 419, 7 

De Grandi, S. & Molendi, S. 2001, ApJ, 551, 153 

De Grandi, S. & Molendi, S. 2002, ApJ, 567, 163 

de Plaa, J., Werner, N., Bleeker, J. A. M., et al. 2007, A&A, 465, 345 



58 



de Plaa, J., Werner, N., Bykov, A. M., et al. 2006, A&A, 452, 397 
Diehl, S. & Statler, T. S. 2006, MNRAS, 368, 497 
Domainko, W., Mair, M, Kapferer, W., et al. 2006, A&A, 452, 795 
Donahue, M., Voit, G. M., Scharf, C. A., et al. 1999, ApJ, 527, 525 
Dressier, A. 1980, ApJ, 236, 351 

Dressier, A., Oemler, A. J., Couch, W. J., et al. 1997, ApJ, 490, 577 
Dupke, R. A. & Arnaud, K. A. 2001, ApJ, 548, 141 
Dupke, R. A. & White, III, R. E. 2000, ApJ, 528, 139 

Durret, K, Kaastra, J. S., Nevalainen, J., Ohashi, T., & Werner, N. 2008, Space Science 

Reviews, 134, 51 
Durret, K, Lima Neto, G. B., & Forman, W. 2005, A&A, 432, 809 
Eckert, D., Neronov, A., Courvoisier, T. J.-L., & Produit, N. 2007, A&A, 470, 835 
Ettori, S. & Fabian, A. C. 2006, MNRAS, 369, L42 
Fabian, A. C. 1994, ARA&A, 32, 277 
Fabian, A. C. & Nulsen, P. E. J. 1977, MNRAS, 180, 479 
Fabian, A. C. & Pringle, J. E. 1977, MNRAS, 181, 5P 

Fabian, A. C, Sanders, J. S., Taylor, G. B., & Allen, S. W. 2005, MNRAS, 360, L20 

Fabian, A. C, Sanders, J. S., Taylor, G. B., et al. 2006, MNRAS, 366, 417 

Feretti, L., Gioia, I. M., & Giovannini, G., eds. 2002, Astrophysics and Space Science Li- 
brary, Vol. 272, Merging Processes in Galaxy Clusters 

Feretti, L. & Giovannini, G. 1998, in Untangling Coma Berenices: A New Vision of an Old 
Cluster, ed. A. Mazure, F. Casoli, F. Durret, & D. Gerbal, 123 

Feretti, L. & Giovannini, G. 2007, arXiv:astro-ph/0703494 

Finoguenov, A., Briel, U. G., & Henry, J. P. 2003, A&A, 410, 777 

Finoguenov, A., David, L. P., & Ponman, T. J. 2000, ApJ, 544, 188 

Finoguenov, A., Matsushita, K., Bohringer, H., Ikebe, Y., & Arnaud, M. 2002, A&A, 381, 
21 

Forman, W, Jones, C, Churazov, E., et al. 2007, ApJ, 665, 1057 
Forman, W, Nulsen, P., Heinz, S., et al. 2005, ApJ, 635, 894 
Fukazawa, Y., Makishima, K., Tamura, T, et al. 1998, PASJ, 50, 187 
Fukazawa, Y, Ohashi, T, Fabian, A. C, et al. 1994, PASJ, 46, L55 
Fusco-Femiano, R., dal Fiume, D., Feretti, L., et al. 1999, ApJ, 513, L21 
Gaetz, T. J. & Salpeter, E. E. 1983, ApJS, 52, 155 
Gallagher, J. S., Garnavich, P. M., Berlind, P., et al. 2005, ApJ, 634, 210 
Gao, L., Navarro, J. F, Cole, S., et al. 2008, MNRAS, 387, 536 
Gastaldello, F. & Molendi, S. 2004, ApJ, 600, 670 

George, M. R., Fabian, A. C, Sanders, J. S., Young, A. J., & Russell, H. R. 2009, MNRAS, 
395, 657 

Gilfanov, M. R., Sunyaev, R. A., & Churazov, E. M. 1987, Soviet Astronomy Letters, 13, 3 
Gilfanov, M. R. & Syunyaev, R. A. 1984, Soviet Astronomy Letters, 10, 137 
Grevesse, N. & Sauval, A. J. 1998, Space Science Reviews, 85, 161 
Gronenschild, E. H. B. M. & Mewe, R. 1978, A&AS, 32, 283 

Gustafsson, B., Karlsson, T, Olsson, E., Edvardsson, B., & Ryde, N. 1999, A&A, 342, 426 
Hamuy, M., Phillips, M. M., Suntzeff, N. B., et al. 1996, AJ, 112, 2438 
Hatch, N. A., Crawford, C. S., Fabian, A. C, & Johnstone, R. M. 2005, MNRAS, 358, 765 
Hattori, M., Kneib, J., & Makino, N. 1999, Progress of Theoretical Physics Supplement, 
133, 1 

Henry, J. P. 2004, ApJ, 609, 603 

Henry, J. P. & Briel, U. G. 1995, ApJ, 443, L9 



59 



Henry, J. P., Evrard, A. E., Hoekstra, H., Babul, A., & Mahdavi, A. 2009, ApJ, 691, 1307 

Ikebe, Y., Makishima, K., Ezawa, H., et al. 1997, ApJ, 481, 660 

Inoue, S., Aharonian, F. A., & Sugiyama, N. 2005, ApJ, 628, L9 

Ivanov, V. D., Hamuy, M., & Pinto, P. A. 2000, ApJ, 542, 588 

Iwamoto, K., Brachwitz, E, Nomoto, K., et al. 1999, ApJS, 125, 439 

Kaastra, J. S. 1992, An X-ray Spectral Code for Optically Thin Plasmas (Internal SRON- 

Leiden Report, Updated Version 2.0) 
Kaastra, J. S., Bykov, A., & Werner, N. 2009, A&A, in prep. 
Kaastra, J. S. & Mewe, R. 1994, Journal of the HEASARC at NASA GSFC, 3, 16 
Kaastra, J. S., Mewe, R., & Nieuwenhuijzen, H. 1996, in UV and X-ray Spectroscopy of 

Astrophysical and Laboratory Plasmas p.411, K. Yamashita and T. Watanabe. Tokyo : 

Universal Academy Press 
Kaastra, J. S., Tamura, T, Peterson, J. R., et al. 2004, A&A, 413, 415 
Kato, T. 1976, ApJS, 30, 397 

Kolmogorov, A. 1941, Akademiia Nauk SSSR Doklady, 30, 301 

Kotov, O. & Vikhlinin, A. 2006, ApJ, 641, 752 

Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 

Landini, M. & Fossi, B. C. M. 1991, A&AS, 91, 183 

Landini, M. & Monsignori Fossi, B. C. 1970, A&A, 6, 468 

Landini, M. & Monsignori Fossi, B. C. 1990, A&AS, 82, 229 

Liedahl, D. A., Osterheld, A. L., & Goldstein, W. H. 1995, ApJ, 438, LI 15 

Lin, Y.-T, Mohr, J. J., & Stanford, S. A. 2003, ApJ, 591, 749 

Lloyd-Davies, E. J., Ponman, T. J., & Cannon, D. B. 2000, MNRAS, 315, 689 

Lodders, K. 2003, ApJ, 591, 1220 

Lutovinov, A. A., Vikhlinin, A., Churazov, E. M., Revnivtsev, M. G., & Sunyaev, R. A. 

2008, ArXiv e-prints 
Mahdavi, A., Hoekstra, H., Babul, A., & Henry, J. P. 2008, MNRAS, 384, 1567 
Makishima, K. 1999, Astronomische Nachrichten, 320, 163 
Mannucci, F, Delia Valle, M., & Panagia, N. 2006, MNRAS, 370, 773 
Markevitch, M. 2006, in ESA Special Publication, Vol. 604, The X-ray Universe 2005, ed. 

A. Wilson, 723-727 

Markevitch, M., Forman, W. R., Sarazin, C. L., & Vikhlinin, A. 1998, ApJ, 503, 77 
Markevitch, M., Gonzalez, A. H., David, L., et al. 2002, ApJ, 567, L27 
Markevitch, M., Govoni, F, Brunetti, G, & Jerius, D. 2005, ApJ, 627, 733 
Markevitch, M. & Vikhlinin, A. 2007, Phys. Rep., 443, 1 
Masai, K. 1984, Astrophys. Space Sci., 98, 367 
Masai, K. & Kato, T. 1987, Physics Letters A, 123, 405 
Mastropietro, C. & Burkert, A. 2008, MNRAS, 389, 967 
Mathiesen, B. F. & Evrard, A. E. 2001, ApJ, 546, 100 

Matsushita, K, Belsole, E., Finoguenov, A., & Bohringer, H. 2002, A&A, 386, 77 
Matsushita, K, Bohringer, H., Takahashi, I., & Ikebe, Y. 2007a, A&A, 462, 953 
Matsushita, K, Finoguenov, A., & Bohringer, H. 2003, A&A, 401, 443 
Matsushita, K, Fukazawa, Y, Hughes, J. P., et al. 2007b, PASJ, 59, 327 
Maughan, B. J., Jones, C, Forman, W., & Van Speybroeck, L. 2008, ApJS, 174, 117 
Mazzali, P. A., Ropke, F. K, Benetti, S., & Hillebrandt, W. 2007, Science, 315, 825 
Mazzotta, P., Mazzitelli, G, Colafrancesco, S., & Vittorio, N. 1998, A&AS, 133, 403 
Mazzotta, P., Rasia, E., Moscardini, L., & Tormen, G. 2004, MNRAS, 354, 10 
McCarthy, I. G., Babul, A., Bower, R. G, & Balogh, M. L. 2008, MNRAS, 386, 1309 
Mei, S., Holden, B. P., Blakeslee, J. P., et al. 2006, ApJ, 644, 759 



60 



Mewe, R. 1990, in NATO ASIC Proc. 300: Physical Processes in Hot Cosmic Plasmas, 
39-65 

Mewe, R. 1991, A&AR, 3, 127 

Mewe, R. & Gronenschild, E. H. B. M. 1981, A&AS, 45, 11 

Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. 1985, A&AS, 62, 197 

Mewe, R., Kaastra, J. S., & Liedahl, D. A. 1995, Legacy, 6, 16 

Mewe, R., Lemen, J. R., & van den Oord, G. H. J. 1986, A&AS, 65, 511 

Meynet, G. & Maeder, A. 2002, A&A, 381, L25 

Million, E. T. & Allen, S. W. 2008, ArXiv e-prints 

Mitchell, R. J., Culhane, J. L., Davison, P. J. N., & Ives, J. C. 1976, MNRAS, 175, 29P 

Molendi, S. & Gastaldello, F. 2001, A&A, 375, L14 

Moore, B., Ghigna, S., Governato, E, et al. 1999, ApJ, 524, L19 

Mushotzky, R., Loewenstein, M., Arnaud, K. A., et al. 1996, ApJ, 466, 686 

Mushotzky, R. F. 1984, Physica Scripta Volume T, 7, 157 

Mushotzky, R. F. & Loewenstein, M. 1997, ApJ, 481, L63 

Mushotzky, R. E, Serlemitsos, P. J., Boldt, E. A., Holt, S. S., & Smith, B. W. 1978, ApJ, 
225, 21 

Nakazawa, K., Sarazin, C. L., Kawaharada, M., et al. 2009, PASJ, 61, 339 
Navarro, J. E, Frenk, C. S., & White, S. D. M. 1995, MNRAS, 275, 720 
Nomoto, K., Tominaga, N., Umeda, H, Kobayashi, C, & Maeda, K. 2006, Nuclear Physics 
A, 777, 424 

Nulsen, P. E. J., David, L. P., McNamara, B. R., et al. 2002, ApJ, 568, 163 

Nulsen, P. E. J., McNamara, B. R., Wise, M. W, & David, L. P. 2005, ApJ, 628, 629 

Oboukhov, A. 1941, Dokl. Akad. Nauk SSSR, 32, 22 

Owen, F. N., Eilek, J. A., & Kassim, N. E. 2000, ApJ, 543, 611 

Paerels, E, Kaastra, J., Ohashi, T., et al. 2008, Space Science Reviews, 134, 405 

Peng, F. & Nagai, D. 2009, ApJ, 693, 839 

Peres, C. B., Fabian, A. C, Edge, A. C, et al. 1998, MNRAS, 298, 416 
Peterson, J. R., Kahn, S. M., Paerels, F. B. S., et al. 2003, ApJ, 590, 207 
Peterson, J. R., Paerels, F. B. S., Kaastra, J. S., et al. 2001, A&A, 365, L104 
Petrosian, V. 2001, ApJ, 557, 560 

Petrosian, V. & Bykov, A. M. 2008, Space Science Reviews, 134, 207 

Poggianti, B. M., Smail, I., Dressier, A., et al. 1999, ApJ, 518, 576 

Pointecouteau, E., Arnaud, M., & Pratt, G. W. 2005, A&A, 435, 1 

Ponman, T. J., Cannon, D. B., & Navarro, J. F. 1999, Nature, 397, 135 

Porquet, D., Arnaud, M., & Decourchelle, A. 2001, A&A, 373, 1110 

Pratt, G. W. 2006, astro-ph/0607199 

Pratt, G. W. & Arnaud, M. 2002, A&A, 394, 375 

Pratt, G. W., Arnaud, M., & Pointecouteau, E. 2006, A&A, 446, 429 

Pratt, G. W., Bohringer, H, Croston, J. H, et al. 2007, A&A, 461, 71 

Pratt, G. W., Croston, J. H, Arnaud, M., & Boehringer, H. 2008, arXiv: 0809.3784 

Rasia, E., Ettori, S., Moscardini, L., et al. 2006, MNRAS, 369, 2013 

Rasia, E., Mazzotta, P., Borgani, S., et al. 2005, ApJ, 618, LI 

Rasia, E., Mazzotta, P., Bourdin, H., et al. 2008, ApJ, 674, 728 

Raymond, J. C. 1988, in NATO ASIC Proc. 249: Hot Thin Plasmas in Astrophysics, ed. 
R. Pallavicini, 3 

Raymond, J. C. & Brickhouse, N. S. 1996, Astrophys. Space Sci., 237, 321 
Raymond, J. C. & Smith, B. W. 1977, ApJS, 35, 419 

Rebusco, P., Churazov, E., Bohringer, H., & Forman, W. 2005, MNRAS, 359, 1041 



61 



Rebusco, P., Churazov, E., Bohringer, H., & Forman, W. 2006, MNRAS, 372, 1840 

Reiprich, T. H., Hudson, D. S., Zhang, Y. ., et al. 2008, ArXiv e-prints 

Reiprich, T. H., Sarazin, C. L., Kempner, J. C, & Tittley, E. 2004, ApJ, 608, 179 

Renzini, A. 1997, ApJ, 488, 35 

Rephaeli, Y. 1977, ApJ, 212, 608 

Rephaeli, Y. & Gruber, D. 2002, ApJ, 579, 587 

Rephaeli, Y, Nevalainen, J., Ohashi, T., & Bykov, A. M. 2008, Space Science Reviews, 134, 
71 

Ricker, P. M. & Sarazin, C. L. 2001, ApJ, 561, 621 

Rizza, E., Burns, J. O., Ledlow, M. J., et al. 1998, MNRAS, 301, 328 

Sanders, J. S. 2006, MNRAS, 371, 829 

Sanders, J. S. & Fabian, A. C. 2002, MNRAS, 331, 273 

Sanders, J. S. & Fabian, A. C. 2006, MNRAS, 371, 1483 

Sanders, J. S. & Fabian, A. C. 2007, MNRAS, 381, 1381 

Sanders, J. S., Fabian, A. C, Allen, S. W., et al. 2008, MNRAS, 385, 1186 

Sanders, J. S., Fabian, A. C, Allen, S. W., & Schmidt, R. W. 2004, MNRAS, 349, 952 

Sanderson, A. J. R., Ponman, T. J., Finoguenov, A., Lloyd-Davies, E. J., & Markevitch, M. 

2003, MNRAS, 340, 989 
Sato, K., Matsushita, K., Ishisaki, Y, et al. 2008, PASJ, 60, 333 
Sato, K., Matsushita, K., Ishisaki, Y, et al. 2009, PASJ, 61, 353 
Sato, K., Yamasaki, N. Y, Ishida, M., et al. 2007, PASJ, 59, 299 
Schmutzler, T. & Tscharnuter, W. M. 1993, A&A, 273, 318 
Schuecker, P., Bohringer, H., Collins, C. A., & Guzzo, L. 2003a, A&A, 398, 867 
Schuecker, P., Caldwell, R. R., Bohringer, H., et al. 2003b, A&A, 402, 53 
Schuecker, P., Finoguenov, A., Miniati, F, Bohringer, H., & Briel, U. G. 2004, A&A, 426, 

387 

Schwartz, S. J., Thomsen, M. F, Bame, S. J., & Stansberry, J. 1988, J. Geophys. Res., 93, 
12923 

Serlemitsos, P. J., Smith, B. W., Boldt, E. A., Holt, S. S., & Swank, J. H. 1977, ApJ, 211, 
L63 

Shapiro, P. R. & Moore, R. T. 1976, ApJ, 207, 460 
Shapiro, P. R. & Moore, R. T. 1977, ApJ, 217, 621 

Simionescu, A., Bohringer, H., Briiggen, M., & Finoguenov, A. 2007, A&A, 465, 749 
Simionescu, A., Roediger, E., Nulsen, P. E. J., et al. 2008a, ArXiv e-prints 
Simionescu, A., Werner, N., Bohringer, H., et al. 2009, A&A, 493, 409 
Simionescu, A., Werner, N., Finoguenov, A., Bohringer, H, & Briiggen, M. 2008b, A&A, 
482, 97 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C. 2001, ApJ, 556, L91 
Sparks, W. B., Donahue, M., Jordan, A., Ferrarese, L., & Cote, P. 2004, ApJ, 607, 294 
Spitzer, L. J. & Greenstein, J. L. 1951, ApJ, 114, 407 
Springel, V. & Farrar, G. R. 2007, MNRAS, 380, 911 

Sunyaev, R. A., Norman, M. L., & Bryan, G. L. 2003, Astronomy Letters, 29, 783 
Sutherland, R. S. & Dopita, M. A. 1993, ApJS, 88, 253 

Takahashi, T., Kelley, R., Mitsuda, K., et al. 2009, in Astrophysics with All-Sky X-Ray 
Observations, Proceedings of the RIKEN Symposium, held 10-12 June, 2008. RIKEN, 
and JAXA Suzuki Umetaro Hall, RIKEN Wako, Saitama, Japan., p.356, ed. N. Kawai, 
T. Mihara, M. Kohama, & M. Suzuki, 356-+ 

Tamura, T., Bleeker, J. A. M., Kaastra, J. S., Ferrigno, C, & Molendi, S. 2001, A&A, 379, 
107 



62 



Tamura, T., Kaastra, J. S., den Herder, J. W. A., Bleeker, J. A. M., & Peterson, J. R. 2004, 

A&A, 420, 135 
Tokoi, K., Sato, K., Ishisaki, Y., et al. 2008, PASJ, 60, 317 
Tozzi, P. & Norman, C. 2001, ApJ, 546, 63 

Tsujimoto, T., Nomoto, K., Yoshii, Y., et al. 1995, MNRAS, 277, 945 
Tucker, W. H. & Koren, M. 1971, ApJ, 168, 283 
Vikhlinin, A. 2006, ApJ, 640, 710 

Vikhlinin, A., Kravtsov, A., Forman, W., et al. 2006, ApJ, 640, 691 
Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009, ApJ, 692, 1060 
Vikhlinin, A., Markevitch, M., Murray, S. S., et al. 2005, ApJ, 628, 655 
Vikhlinin, A., Voevodkin, A., Mullis, C. R., et al. 2003, ApJ, 590, 15 
Voigt, L. M. & Fabian, A. C. 2006, MNRAS, 368, 518 
Voit, G. M. 2005, Reviews of Modern Physics, 77, 207 

Voit, G. M., Balogh, M. L., Bower, R. G., Lacey, C. G., & Bryan, G. L. 2003, ApJ, 593, 272 

Werner, N., Bohringer, H., Kaastra, J. S., et al. 2006a, A&A, 459, 353 

Werner, N., de Plaa, J., Kaastra, J. S., et al. 2006b, A&A, 449, 475 

Werner, N., Finoguenov, A., Kaastra, J. S., et al. 2008, A&A, 482, L29 

Werner, N., Zhuravleva, I., Churazov, E., et al. 2009, MNRAS submitted 

White, D. A. 2000, MNRAS, 312, 663 

White, S. D. M., Briel, U. G., & Henry, J. P. 1993, MNRAS, 261, L8 
Wik, D. R., Sarazin, C. L., Finoguenov, A., et al. 2009, ApJ, 696, 1700 
Woosley, S. E., Kasen, D., Blinnikov, S., & Sorokina, E. 2007, ApJ, 662, 487 
Xu, H., Kahn, S. M., Peterson, J. R., et al. 2002, ApJ, 579, 600 
Zhang, Y.-Y, Finoguenov, A., Bohringer, H., et al. 2007, A&A, 467, 437 
Zwicky, F. 1937, ApJ, 86, 217 



